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O ■ Abstract 



We consider the quantum dynamics of many bosons systems in the mean field hmit with 
a singular pair-interaction potential, including the attractive or repulsive Coulombic case in 
' three dimensions. By using a measure transportation technique developed in [3], we show that 

Wigner measures propagate along the nonlinear Hartree flow. Such property was previously 
proved only for bounded potentials in our works [51 [B] with a slightly different strategy. 
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1 Introduction 



The evolution of a system of many quantum non relativistic particles is described by an n-body 
Schrodinger Hamiltonian. The mean field limit consists in replacing this problem by a non linear 
\ 1-particle problem, by considering that one generic particle interacts with the average field of all 

particles, in the limit when the number of particles is large and the interaction potential is inversely 
\ weak. Practically, it is common knowledge that this approximation starts to be very effective when 

• the number of particles exceeds a few tens. In the last decades, many works have been devoted 

to justify this limit. Most of them considered the mean field dynamics of well prepared quantum 
states, coherent states or Hermite states, by following and extending the phase-space approach, 
also known as the Hepp method (see I151[171[311[33IM1I11[S7]) or by studying the BBGKY 
hierarchy of reduced density matrices (see [H [171 [HI ISOl HH [59]). Some of these results achieved 
to include very singular pair interaction potentials in [SI [HI HOI HI] or considered the rate of 
convergence (see [7] [571 HI] ) , sometimes motivated by the modelling of Bose-Einstein condensates 
. , (see a.e. |21[ 1^51 [Tj). In this article, we continue our program, which consists in deriving the mean 

' field limit, for general initial data in the bosonic framework. Our strategy is inspired by older 

attempts to give substance to the formal link between bosonic Quantum Field Theory and the 
finite dimensional microlocal or phase-space analysis (see [ini [IS HH HSl HZ] ) • With this respect, 
the small parameter £ = asymptotics is the infinite dimensional version of semiclassical analysis. 
And it has been realized in the 90's, that the Wigner (or semiclassical) measures provide a powerful 
tool in order to obtain the leading term in the semiclassical limit (see [291 [5D1 1571 HI]), because 
they flexibly and efficiently incorporate a priori estimates (see [T51 [HI [^ [511 153| ). 

In [4] Wigner measures were introduced in the infinite dimensional setting and their main 
properties were studied. This presentation exploited and clarified the intimate relationship be- 
tween pseudo-differential calculus, phase-space geometry and the probability approach, inherent 
to bosonic QFT. In 5 , the dynamics for well prepared data and bounded interaction potentials 
was reconsidered within this approach. The general propagation result was obtained in ^ for 
bounded interaction potentials. In particular, we showed that the BBGKY hierarchy dynamics is 
a projected picture of the evolution of the Wigner measure, for which there is a closed equation. 
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One difficulty which was solved in [5^ is concerned with the integration of a weak Liouville equa- 
tion valid after testing with cylindrical or polynomial observables: Such classes of observables are 
not preserved by the nonlinear Hamiltonian mean field flow. For bounded interaction potentials, 
the number conservation allows polynomial approximations of the nonlinear deformation in balls 
of the phase-space. This done by following a truncated Dyson expansion approach presented in 
[251 Ell '^27^ In applications, an important case is the one of 2-body Coulombic interaction because 
it models the general non relativistic motion of charged (or gravitational) particles. Again there 
are results about the mean field problem for specific initial data (see [SI [44]), but the approach 
that we followed in [5] essentially fails. With a singular pair interaction potential, a solution to 
this problem is provided by measure transportation techniques developed for optimal transport 
theory (see [31ini]). Hence, the dynamical mean field limit relies even more on the fact that Wigner 
measures are probability measures on the phase-space. 

We now expose our main result. The Hamiltonian of an n-body quantum system, with a pair 
interaction potential, is given by the Schrodinger operator 

n 

i—l l<i<j<n 

where e is a positive parameter and Xi,Xj € R'^. We assume that the particles obey Bose statistics. 
So, we consider ifi"-* as an operator acting on the space L^(R''") of symmetric square integrable 
functions. This means that 

*eL2(]^dn) iff ^' e i2(-jjdn^ ^'(xi,--- ,x„) = *(a;^i,...,a;^J a.e 

for any permutation a on the symmetric group ©„. The mean field asymptotics is concerned with 
the limit as e — >■ and — >■ 1, where represents the number of particles of the system. 

Let H be the direct sum of Hilbert spaces of the form 



H = 0i?(M'i»), 

n=0 

and consider the Hamiltonian of the many-bosons system (with arbitrary number of particles) as 

C30 

if,=0i?i"). (1) 

ji=0 

An obvious feature of the operator is the conservation of the number of particles. Hence, it is 
useful to define the number operator 

oo 

N = 0£n ll2,2(Rd„). 

n=0 

The free Hamiltonian, corresponding to = 0, will be denoted by H^: 

oo n 
n=0 i=l 

Second quantization is a natural framework for the study of many-body problems and, even more, 
it helps to understand the mean field limit and the structures behind it. However, the result can be 
presented without using the language of quantum field theory. We just mention that the operator 
Hg can be formally rewritten as 

He— W a* (x) .W a{x) dx + - / V{x ~ y)a* {x)a* {y)a{x)a{y) dxdy , 

.iBd 2 .lo2d 



2 



with the £-dependent canonical commutation relations [a{x) , a*{y)] — sS{x — y) ■ It is interpreted 
as the Wick quantization of the classical Hamiltonian 

h{z,z)= f \Vz{x)\^ dx + \ [ \z{x)\Mv)?V{x-y) dxdy. (2) 

In our analysis an operator, which violates the number of particles conservation, will play an 
important role, namely the Weyl operator. Such operators are given for / G Lp'iW^) by 

W{f) = eTf , 

where a*{f),a{f) are the creation-annihilation operators on H satisfying the e-canonical commu- 
tation relations (CCR): 

Hh),a*{h)] = e{h,h)mM.-) X [a*(/i),a*(/2)] = = [a(/i), a(/2)] . 
Accurate definitions on second quantized operators can be found in Appendix [Bl 

Our approach is based on Wigner measures which are Borel probability measures on the infinite 
dimensional phase-space Zg := L^(R'',C). The states of the many-bosons system are positive 
trace-class operators on % of normalized trace equal to 1 (i.e.: normal states or density operators). 
To every family of those states {Qe)e&{Q,e) we asymptotically assign, when e — > 0, at least one 
Borel probability measure /i on Zq := L^(K'^,C), called Wigner measure, such that there exists a 
sequence (efc)feeN, limfe E/t = and 

hm TrK, W{V2i:0] = , 

A:— >u 

under the sole uniform estimate Tr [fPeN*] < Cs for some 5 > 0. Here !F~^{^) is the inverse 
Fourier transform of /i. 

The problem of the mean field dynamics questions whether the asymptotic quantities, namely 
Wigner measures, as e — > associated with 

Q^{t) ^ e-'^-"- Q^e^i"' , isM 

are transported by the flow s) generated by the classical Hamiltonian h{z, z) and given, after 
writing Zt = ^{t,s){zs), by 

idtzt = {d-,h){zt,zt) = -Azt + V * \zt\''zt . (3) 

After checking that the Hamiltonian ^ has a self-adjoint realization so that the quantum dynamic 
is well defined on T-L and after checking that the mean field flow is well defined on Zi ~ H^{W^), 
our main result is stated below. 

Throughout the paper we assume that the real valued potential V satisfies the assumptions 

V{-x) ^V{x) <=,R, (Al) 
V(i - A)-i/2 e C{Za) , (A2) 
and (1- A)-i/2-^(l- A)-i/2 G £°°(Zo). (A3) 

Here we use the notation £(()) for the space of bounded operators on the Hilbert space f) and 
'C^(f)), 1 < p < +00, for the Schatten classes, C°°{\)) being the space of compact operators for 

P — +CO . 

Theorem 1.1. Let {gi;)^^(^Q g-^ be a family of normal states on "H with a single Wigner measure hq 
such that the bound 

Tr[(N + i/O) Ve] < < +0O , (4) 
holds uniformly w.r.t e G (0, e) for some (5 > . 

Then for all t G M, the family (e~*~^=£i£e*~^=)£g(o,£) ^'^■s o, unique Wigner measure m which is a 
Borel measure on Zi = H^{K^) . This measure fif ~ $(t, 0)*/io is the push forward of the initial 
measure /io by the flow associated with well defined on Z\ . 
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In a formal level the proof of the latter theorem is rather simple. We first write an integral 
formula 

ft 4 

Tv[QS)Wm=Tv[e,Wm+i / Tr[^.,(5)W^(e)E^'"'^^]^^' 

where Qe{t) = e"**/*^^^ gi^e**/^^^ and Oj are some Wick quantized observables. Then we attempt 
to take the limit as e — , where in the r.h.s only the term with j = 1 will contribute. We end up 
in a formal level with a transport equation on the Wigner measures 

dt^i + i{h , ^} = Q, {h, fi} = dzh dgfJ, — dzH dgh 

which we then solve by appealing to the results in [3]. 

Outline: In Section [2] we prove the self-adjointness of the Hamiltonian and the existence of 
a global flow on Zi = H^{M.'^) for the Hartree equation The derivation of the mean field 
dynamics is done in Section [3] where we prove Theorem ll.il Some additional properties are stated 
in Section |4j in particular, we draw the link with former results on bounded potential and reduced 
density matrices and provide non trivial examples elucidated by the Wigner measure approach. We 
finish the article with several appendices dedicated to second quantization, absolutely continuous 
curves in Prob2(2) as well as some weak conditions for the potential V ensuring the fulfillment 
of the assumptions and {J^. 

2 Well defined dynamics 

In this section we shall prove that: 

• the quantum dynamics is well defined, namely iJ^ has a natural self-adjoint realization; 

• the mean field dynamics is well defined on Zi ~ H^{M.'^) , with additional useful estimates. 

2.1 Self-adjoint realization of 

The Hamiltonian has a particular structure explained in a general framework in Appendix |^ 
Let y be a real- valued Lebesgue measurable function a.e. finite satisfying assumptions (A[T]) and 
{J^. The multiplication operator 

l<'i<j<n 

with its natural domain ^(Fi"^) = {* G i^^j^dn) . yin)^ ^ L^R'^'')} is self-adjoint on 
as well as the differential operator 

n 

fjUn) = £^ _A:,, , with V{H°^^"^) = LliR''") n ij2(R'^«) . 

i=l 

Therefore, according to Appendix \X\ 

oo oo 

V,^J2^t\ and i/° = 5: 

n=0 n=0 

endowed with their natural domains are self-adjoint on H. 

Proposition 2.1. Under the assumptions and (J^: 

(i) The operator 
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is self-adjoint on V{He'^"') C V{V}">). 
(ii) The operator 

oo oo 

n=0 n=0 

is self-adjoint and essentially self-adjoint on ®n=o'^n where Dn is any core of Hs"'\ 

Proof, (i) By assumption we see that v}"'^ is infinitesimally small with respect to H^'^^\ 

So that, V{He'''''^) C 2?(K^"^) and the operator ijj"^ = iJ^"^^"' + yj"' is self-adjoint on the domain 
of 77°'^"-* by Kato-Rellich theorem. 

(ii) Applying Proposition lA.il we see that is self-adjoint and essentially self-adjoint on ©°'f q^"- 
□ 



Later, it will be useful to use the reference operator 

oo 

S,{X) = ^ + en + X{enf . (5) 



n=0 

which is self-adjoint by Proposition lA.ll Moreover, by functional calculus of strongly commuting 
self-adjoint operators we observe that 2?(S'e(A)) is invariant with respect to the parameter A > 0. 

Proposition 2.2. Under the assumptions and (J^, for any A > 0, the operator is 

Se{X) -bounded with 

v* e v{s,{X)), WVM-H < M\v{i - A)-^/^cimm) \\Se{x-^)n-H ■ 

Therefore is essentially self-adjoint on T>{Se{X)) . 

Proof. The multiplication operator by V{xi — X2), at least defined as a symmetric operator from 
into 5'(]R2d)^ satisfies 

- T/(xi)(l-A,J-i/2e/:(L2(K2rf)). (6) 

For ^ G "H, $ G P(5e(A)), taking advantage of the symmetry of those wave functions, we compute 

00 

{^,V,^) - ^(*("),n(n-l)£V(a;i-X2)$("))^.(R.„) 

n=2 
00 

= , n{n - l)eV(xi - X2){1 - A,J-V2(i _ A, Ji/2$(«))^,^^,„^ . 

n=2 

By noticing that 

n 

1=1 

= ||N-''/2(N + i/0)l/2$(«)||2^, 

the Cauchy-Schwarz inequality leads to 

KM/, K<f>| < |ini-A)-i/2|| ||vi,||^ WN'/^N + H^)'/H\\n ■ (7) 
Now with the inequality ab < (Aa)^ + (&/A)^, we see that 

||N3/2(N-f ff^")i/2$||2^ = ($,N3(N-KffO)$) 

< ($,A-2N^ + A2(N + iJ^")2$) 

< A2||5,(A-2)$||?,. (8) 
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Putting together ([7]) and ([5]) yields the estimate. □ 

We end this section with some invariance properties of the domain ^(^^(A)) with respect to 
the Hamihonian and the Weyl operators. 

Proposition 2.3. For any A > and t eM. 

e-»7^^I?(5,(A)) cI?(5e(A)). 

Moreover there exists Cx > such that 

\\S,{X)e-'i"^{S,{X) + ir^Wan) < Cx for all t e R. 

Proof. For any e ^(^^(A)) C V{H^), we observe that e~*^^-* belongs to V{N^)nV{Hs) since 
I?(5e(A)) C V(N^) and strongly commute with N. Proposition O implies + K) = 

ViSsllS)) = T>{SeiX)) when /3 > is large enough, so that for any $ e V{Se{X)) 

and hence e-*t^-* belongs to V{{Se{P) + V^)*) = V{SeiP)) = T^iSeW)- 
Again for /3 large enough, we know that 

1 + SM + v,^{i + v,{SM + ir^){s,{p) + 1) , 

and 

(1 + N + /3N3 + H,)-' = [SM + + V,{S,{P) + 1)-')-' ■ 

Therefore the operators (1 + + N + I3N^ + H^)-^ and {1 + N + I3N^ + H,){SeW) + l)-^ 

are bounded. We conclude with 

5,(A)e-'i^' (1 + S.iX))-' = S,{X){1 + SMr'i-^ + + N + /JN^ + H,)-' 

o e-'i"^ (1 + N + /3N3 + H,)il + S,{(3))-\l + S,il3)){l + 5e(A))~' . 

□ 



Proposition 2.4. For any ^ e H'^{W^) and A > the domain 2?(S'e(A)) is invariant with respect 
to Weyl operators W{^) with 

\\{S,iX) + ir'W{OS,{X)\\cin)<Cx,i, 
uniformly w.r.t e £ (0,£) for some constant Ca,^ > . 
Proof. For ah $, e X>(5e(A)), we write 

($, W{0*S,iX)Wim = (-f, (5,(A) + QY^'^m, 
where is the fohowing polynomial 

IS IS ts 

Q,{z) - (Z + —i, -A{Z + -^0)20 - -Mzo + Pe{z + ^0 - Pe{z) 

and Peiz) — + 3£|z|2(, + s^NI^q is the complete Wick symbol of N'^, according to Propo- 

sition [B]2] or by direct computation. The assumption ^ e ensures that is uniformly 

bounded in (Sp+qK^Vp.qiZQ) and the number estimate of Proposition IB. 31 says that QY^'^'^{'^)~^ 
is a bounded operator and therefore 

Hence for * e ^(^^(A)), 



I nWick 



iSeiX) + 



(5e(A) + l)^' 



and VK(0* belongs to V{Se{X)) , with \\Se{X)W{0'f\\ < CxM^eW + l)*ll ■ 1^ 
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Proposition 2.5. For any function x G C^(R^) and A > 0, the operator x(N, H^) satisfies 

Vfc e N, ||N^-5aA)x(N,i?,)||£(„) < Cl+^ 

for some C\^^ > . 

Proof. The operators N, (like N and (N + H^) are strongly commuting self-adjoint operators 
so that the functional calculus is well defined for the pair (N,_ffj) . With a cut-off function xi G 
C^(K.) such that Xi(cc) = 1 on a neighborhood of supp x, the operator N*''(l + N + /SN'' + 
iJ£)xi(N)x(-ffe,N) is bounded with 

II (1 + N + /3N3 + i?,)N'=xi(N)x(i/, m\cm < CpCl+^ . 

For sufficiently large /?, Proposition 12.21 savs 

\\S,mi + N + + H.r^Wc^u) < C'p . 

This is done with 

N'=5,(A)x(N, H,) - S,{\){1 + + 5e(/3)(l + N + /SN^ + i/e)"'x(N, i/^) , 

andCA,x =niax{Q, CpC'^\\S,{X){l + S,{(3))-^\\Y □ 



2.1.1 Mean field dynamics 

We shall use another more convenient writing of the Cauchy problem 

idtzt = -Azt + V * IzA'^Zt 



(9) 

zt=o = zo- 



After setting zt — e**' ^''zt — e ^*'^zt it becomes 

idtit = e-'*^ [V *\e''^zt\'^{e'*^zt)\ 



2:4=0 — Zq . 

Proposition 2.6. Assume o,nd (j^. For any zq E Zi = H^{M.'^) the Cauchy problem ^ 

admits a unique solution {t zt) S C°(R; i/i(IR'')) nCi(M; H-'^{K^)) . More precisely, the Cauchy 
problem (jlOp which is equivalent to ^ admits a unique solution in {W\ {W^)) . Moreover these 
solutions verify 

\Zt\L^ = \zt\L^ = \zo\l^ (11) 

and h{zt,z[)^h{zQ,z^), (12) 
for h{z,z)^ f \VzWx)dx + \[ V{x-y)\z{x)\Mv)? dxdy. 

Finally, the time- dependent velocity field defined on R x Zi by 

vit,z) = e-''^i[V*\e''^z\^]e''^z) 

satisfies the estimates 

\vit,z)U„ <\\V{1- Ar'/'\\\z\%JzU, (13) 
and \vit,z)\z,<\\V{l-A)-'/'\\\z\l^\z\z„. (14) 

Proof. The first results are standard (see a.e. [121 131]) in the analysis on nonlinear evolution 
equation. Nevertheless, we recall the details of the proof be cause it also contains (|T3)) (IT4l) . which 
is crucial in our analysis. 



7 



By considering the second formulation pUj) . it suffices to prove that the mapping z —i' (V * \z\'^)z 
is locally Lipschitz in H^{M.'^) . After noticing that the distributional derivative of {V * \z\'^)z is 

dx[{V *z^Z2)z3] = V* {dxZiZ2 + zidxZ2)z3 + {V *ziz2){dxZ3) , (15) 

it is reduced to the estimate oi V * (2122)23 in in terms of the and H^- norms of 21, 22, 23 . 
For C € i^(K''), write 

{^,V* (2122)23) i2(Rd) = (Zi (g) ^ , V{xi - X2)Z2 ® 23) L2(R2d) . 

When b is the multiplication operator by V{xi — 2:2), the estimate (O says that 6(1 — A^j^)"^/^ is 
bounded, with 

\{^®zi,V{xi-X2)z2®Z3)l2(s,2.)\ < ||1/(1- A)-i/2|| \S_\^2\zi\l2\z2\hAz3\l^ ■ 
A symmetric argument of fG]) says b{\ — A^jj)^^/^ is bounded, with 

1(^(^21, ^(Xl - 0:2)22 (»23)i,2(R2d) I < ||V^(1 - A)-l/2|| |^|^2|2l|i2|22|LH23|ifl • 

Finally the symmetry of the expression V * (2122)23 w.r.t exchanging zT and 22 gives 

1(^0 21, F(a;i-X2)22® 23)i2(R2.) I < |ly(l- A)-i/2|| |eliH^l|i/ik2|LH^3|L2 ■ 
We have proved, owing to (fTS)) . 

(2122)23^2 < ||V^(1- A)-i/2|| mill |2,(i)|Hi|2,(2)|L2|2,(3)|i2 (16) 

from which we deduce 

|V^* (2122)231^1 < 11^(1 -A)-l/2|, ^^^^ |2,(i)|ffl|2,(2)|ffl|2,(3)|L2. (17) 

Since z — e'*^z preserve the and i/^ norms, the velocity field estimates and (|14l) are 
consequences of (ITB|) and pT|) . 

For the sake of completeness, let us finish the proof of the global well-posedness of the Cauchy prob- 
lem. The estimate (fT7|) provides the Lipschitz property of 2 — >■ * |2p2 in H^iW^) . This implies 
the local in time existence and uniqueness of a solution to (fTO|) in C"'^((— Tz^, Tz^); 77^(R'')) , and 
therefore the local in time existence and uniqueness of a solution to ([9|) C°((— Tz^, Tj^); H^{W^)) n 
Ci([-rzo,rzo];i?"i(M'')) . The global in time existence then comes as usual from the control of 
12*1^1 — \zt\m deduced from the conservations of ([TT|) and (fT2|) . For ((TT|) . take the real part of 
the scalar product of each member of (O with It . This implies 9t|zt|^2 = . 
For (|12p take the scalar product with x(— i?^^ A)9t2t where x G C^(R) satisfies < x ^ 1 and 
X = 1 in a neighborhood of 0, with i? > : 

= 2Re {dtzt , x(-ii'"'A)i9t2t) 

= 9t(2( , -Ax(-i?-^A)2t) + 2Re {dtZt , x(-i?"'A)[C^' * \zt?)zt]) 

Integrating this identity from to t and taking the limit as, co with the help of (fT7| gives 
f |V2t|2dx - / IV20P da; + 2 / Re (9,z, , {V * |z.f )2,) = . 

JWL'^ JR'' Jo 

Due to the symmetry of V{x) = V{—x), the last integrand equals 

Re {dsz, , (V * 12,12)2,) - / a,(|2,(a:)n\/(a; - y)\zs{y)\^ dxdy 

= / \zs{x)\^Vix-y)\z4y)\^ dxdy. 

The conserved quantities ([TT|) and (IT^ combined with imply |zt|jji < C|2o|//i for some con- 
stant independent of t G (— T^q, T^^,), and hence Tz„ = +00 . □ 
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3 Derivation of the mean field dynamics 

This section contains the proof of our main Theorem ll.il Below, we recall from our previous work 
[3] the notion of infinite dimensional Wigner measures and collect some of their properties. We 
will often make use of Weyl and Wick quantization throughout this section. So, we suggest first 
the reading of Appendix [Bl 

Two phase-spaces will be necessary for this analysis: Zq = L'^{W'';C) (resp. Zi = iJ^(R'';C)) 
endowed with its scalar product ( , ) (resp. {zi,Z2)zi = {zi , (1 — A)z2)), its norm \z\'^^ = 
{z , z) = jzlla (resp. = NllfOj its real scalar product Re (z , z) (resp. Re (zi , Z2)zi)- Only 

on Zq, we will use the symplectic structure with a{zi,Z2) = Im (zi , Z2) . Meanwhile, the real 
euclidean structure on Zi is important especially when referring to Appendix [Cl 

3.1 Wigner measures 

The Wigner measures are defined after the next result proved in [4] Theorem 6.2]. 

Theorem 3.1. Let {ge)^^(^Qg^ be a family of normal states on % parametrized by e. Assume 
TtIq^'N^] < Cs uniformly w.r.t. e € (0,e) for some fixed 6 > and Cg G (0, +00). Then for every 
sequence (e„)„gN with lim„_^oo £n = 0, there exist a subsequence (£„^)fcgpf and a B or el probability 
measure on Zq , such that 

lim Tr[£-g„ b^'^y^] = 

for all b G Scyii^o). 

Moreover tMs probabMy measure , satisfies f \z\% d,iz) < 00. 
Definition 3.2. 

The set of Wigner measures associated with a family {Qe)ee(a,e) (resp. a sequence {gf;^)nefi) which 
satisfies the assumptions of Theorem \3.1\ is denoted by 

M{ge,e e (0,e)) , {resp. M{ge„,n e N)) . 

Moreover this definition can be extended to any family {ge)ee{Q,e) such that 

||(l + N)Ve(l + N)^||£i(«) <C5 

for some S > with the decomposition g^ — Xf '^gf '^ ~ Xf^'^ gf '^ + iX^'^ g^'^ — iX^'^ gl'~ . 

Wigner measures are in practice identified via their characteristic functions according to the 
relation 

M{ge,ee{0,e))^{fi} ^ \imTT[g,W{V27:0] = T~\pm 

O \imTT[g,WiO]= [ e''^^<= <«'^>dA<(^). 

The expression A4{ge,s £ (0,e)) — {n} simply means that the family {ge)Ee{o.e) is "pure" in the 
sense 

limTr [g.b^^y^] = J b{z) dfi , 

for all cylindrical symbol b without extracting a subsequence. Actually the general case can be 
reduced to this one, after reducing the range of parameters to e € {Sn^ , fc S N} . For checking 
properties of the elements of A4(g^,e G (0,£)), extracting a subsequence in this way allows to 
suppose without loss of generality A^(gie, e G (0, e)) = {/i} . 

A simple a priori estimate argument allows to extend the convergence to symbols which have 
a polynomial growth and to test to Wick quantized symbols with compact kernels belonging to 
^a%(2o) = (St^f.'P^.iZo) (see H Corollary 6.14]). 



/ Hz) d^i{z) , 

JZc, 
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Proposition 3.3. Let {ge)^^^f) ^) be a family of normal states on C(T-L) parametrized by e such 
that Tr[geN"] < Ca holds uniformly with respect to e ^ (0,£), for all a G N,and such that 
M.{Qe,s G (0,e)) — {/i}. Then the convergence 

lim Tr[£-e&^'^''] = / b{z) dpL{z) , (18) 
Jzo 

holds for any b G 'P^giZo)- 

Wigner measures are completely identified via Weyl-quantized observable. One may question 
whether testing on all the ^^^^'^ with b G V^giZo) also identifies the Wigner measures. Indeed, 
this is the case according to the result proved in [4j Theorem 6.13]. 

Proposition 3.4. Assume that the family of operators {Q^)ee(Q.e) satisfies 

||(l + N)>^(l + Nn|^,(^) <C„ 

uniformly w.r.t e G (0,e) for all a G N. For any fixed /3 belonging to 'P'^g{Za) the family 
{13^^'^'^ Q^)eG{o,e) satisfies the assumptions of Definition \3.'2\ and 

M{p'^^'^^g-)^{P^i, ^ieM{Q')} . (19) 

3.2 Weak mean field limit of the dynamics in terms of the characteristic 
function 

After some extraction process and for some specific initial data {Qs)ee(a.E)^ ^ family {iJLt)teWL of 
measures can be defined and solves weakly a transport equation. We consider on L^(R^'') the 
(unbounded) multiplication operators 

V = \v{xi - X2) and Vs = [e-''^^^ ® g-^^A., )y (e'''^==i ® e'^^==2 ) , 

and respectively associate with them the polynomials, well defined on Zi = i/^(R''), 

y(z) = i(z«2^V^(x-y)z«2)i.(K2.) and V;(z) = (z^^^ z^^^^.^j..) , z G . 

Instead of considering 

gS) = e-'i"^ Bec'i"^ . 

we will rather work with 

~g^{t) = e^i"'^e-^i"^e,e^i"^e-^i"° . (20) 

Our assumptions will be made in terms of the operator 5'e(l) already introduced in ([5]) and which 
can be rewritten with Wick observables. 

Definition 3.5. The operator Se is defined by 

oo 

Se = Y^ + + i^nf = dr(i - A) + N3 , 

ra=0 

with domain V{S,) = j^GH, E^^o + + (e")^)**"^ lli2(Rd„) < oo} and H°'^'''^ = 

dr(i-A)lv"^„- 

Remember that it is self-adjoint with this domain (see (O). Moreover it can be written Sf. = 

s,{z) = (z , (1 - A)z) + [|z||„ + 3e|z||„ + sMU ■ 
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Proposition 3.6. Let {Qe)^^{j;) &e a family of normal states on % satisfying for some constant 
C > the estimate 

Tr[(l + Se)Qe{l + Se)] < C uniformly w.r.t e £ (0,e) . 

The operator is the one given in Definition \3.5\ and (t) is the operator given by (|20|) . Then 
for any sequence (£n)nGN in (0,£) such that lini„_j.oo £n = there exist a subsequence (e„^)fegN and 
a family {flt)teR of Borel probability measures on Zq satisfying for any i € M 

X(fe,.,(t),fceN) = {At}, 

with the Liouville equation 

-^(gZV^Ro (c,.)) ^ ~^(g.y2Rc(«,.)) _ 272* /* -^sie'^^' Im (z®^ V;^ ® z)) ds , 

Jo 

= Ao(e^^^° / M.({T4(.);e'^''' ^^'-^Dds, for allS, £ . {21) 

Jo 

Proof. The proof uses several preliminary lemmas which we defer below. The first step is to prove 
the existence of Wigner measures defined for all times t € R . This is done in Proposition [2111 Let 
us now prove the Liouville equation. 
By Lemma 13.81 we have 

Tv[~gS)Wm = Tr[^,,W^(0] + ^ / Tr[&(s)M^(0 V e^-^6,(s, 0"^^^'] ds , (22) 

where bj are the following polynomials 

= -2V2 Im V; C ® z) 62(5, = -Re {z'^^ t C®') + 2(e V z, V;C V z) 

63(5,0 = y2im v; z) 64(5,0 = v; ^^^^ _ 

With the number estimate in Proposition IB.H Lemma 13.71 below will ensure that the sum in the 
r.h.s over j = 2, • • • ,4 converges to when e 0. On the other hand, the term with j = 1 has a 
limit according to Lemma r3 . 1 01 applied with ge{t) after noticing that Tr [(1 + Sf;)ge{t){l + S^)] < C 
owing to 11(1 + Se)^^e'i"oe-'i"-{l + S'e)=Fi|| < C" due to Proposition O □ 

The above proof is completed in essentially three steps: 1) The relation (|22|) is first established 
by extending Wick-calculus arguments to the case when V is unbounded, and rough estimates for 
bj{s, ^^^'"^''^ j — 1, ... ,4, are given; 2) An Ascoli type argument, relying on these rough estimates 
allows to make the subsequence extraction [eni,)keN uniform for all t G R; 3) An additional 
compactness argument is given in order to ensure the convergence of the term with j = 1 in p2p . 

3.2.1 Wick calculus with unbounded kernels 

The results presented in this paragraph would be direct applications of the Wick calculus given in 
Proposition [B2] for a bounded potential V G L°°(R'^) . Although the algebra is the same as in the 
bounded case, justifying the formulas for unbounded potentials fulfilling (AH]) (Al2]) (A|3]) requires 
some analysis. 

Lemma 3.7. The identity 

(V.^'^'^WiO - W{OV,^'^^) * = W{0 ^^e'6j(s,0^"''j * • (23) 

holds for any ^ G H'^(R'^) and ^ € T>{Sg), with given by Definition \3. 51 Additionally, for all 
^' e V{Se) C V{Sl''^) C X>(Ni), the estimates 

\\b,{s,s,r^^^ n < c(i + leilJIKi + N)i*|| < c'ii + i^iijiKi + s,y/H\\ 

<C'(l + |^|ij||(l + 5,)vI/||, (24) 
hold uniformly w.r.t j G {1, . . . , 4}, s G R, when ^ e iJ^(R'*) and therefore when ^ G H^{R''-) . 
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Proof. We first remark that, owing to the assumption and the estimate ([TB)) . the polynomials 
bji^^O: j — ■ ■ A belong to the set ®p^q<3Vp^q{Zo) , with 

l^jlep+,<3-Pp,,(2o) ^ + iClli) • 
Hence, Proposition [B?T] and Proposition 12.21 prove ([M]) with 

I?(5e)Cl?(N3/2)c2?(6j(s,0^"''), J = l,---,4, and C PlV^^^'^) . (25) 

By Proposition HH] the domain I?(S'e) is invariant with respect to W{^) for aU ^ e i/^(M'^) . A 
Taylor expansion yields, for all z e Zi, the equality 

4 

The formula ^ is standard for bounded V due to W^*(0&'^*'^''H/(0 = K + ^^)^"''' when 
^ G 'Paig{2o) = ffip'qgN^p.9('^o) • Let US rccousidcr the proof of this result for our unbounded V . 
With the previous estimates, the quantity A{t) = {'^,W{tC)Vs{. + ^t^)^'='=iy(iO**> is well 
defined for all * e V{Se) with 

4 

We first establish in a weak sense the equality (|23l) : Differentiate for any T^{Se) , 

3 

where ho{z) = Vs{z) . Now, a direct calculation with 0(^) = + gives 

for j = 0, • • • , 3 . Therefore .4(1) = .4(0) and, knowing (P5|) . we conclude that 

W^(0 I y^""^ + X! ^'^J^*' ^)^"'' I = . (26) 

for any e I?(5'e). Now taking = PF(0^ in ([Ml) for any ^ e X'(5e), while * G ^'(S'e) owing to 
^ e i?^(M''), we obtain the claimed equality. 

□ 



Lemma 3.8. Let {Qe)e£{Q,e) family of normal states on H. Assume that gdSe + 1) G L^{H.) for 
all e G (0,e), with given by Definition \3.5\ and Qe{t) by (l20l) . Then for any ^ G H'^(M.'^) , the 
map s i~> Tr[^e(s) VF(C)] belongs to C^(R) and the following integral formula holds true 

^ "'0 3 = 1 
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Proof. We write 



Tr[(^e(s) - ge{s))W{0] 



= Tr 
+Tr 



The following limits hold true on ^{Se) 
1 



lim 



lim 

s^t t — S 



by Stone's theorem and the invariance of 'D{Se) w.r.t e'*'^" and e'*^^ . By using the estimate 
in Proposition 12.31 we see that the latter limits are limits in C{Ti.) w.r.t the strong convergence 
topology. After noticing that g^e^'^^ e~^~^' W{^){S£ + 1) is trace class when ^ g i/^(R) owing to 
Proposition [23] and Proposition 12. 3) we are done by taking the trace and letting s — > i. 
Now it suffices to take the integral to get 

£ Jo 

When ^ e H'^(&.'^) , the equality 

4 

Tr [(1 + SMs) {V^^'^WiO - W^(OK^^^') (1 + S,)-'] ^ TT[g,{s)W{0 E 
makes sense, since (1 + 5e)^e(s) G C^{H) and 

by Lemma 13.71 □ 



Wick^ 



i=i 



3.2.2 Subsequence extraction for all times 

The first step in the proof of Proposition 13.61 is to show the existence of Wigner measures for all 
times. This is accomplished below following merely the same lines as [S!, Proposition 3.3]. 

Proposition 3.9. Let (£ie)^g(-Q be a family of normal states on % satisfying for some constant 
C > the estimate 

Tr[gie(l + Se)] < C uniformly w.r.t e £ (0,e) . 

The operator and ge{t) are respectively given by Definition \3.5\ and (j20p . Then for any sequence 
(£n)neN in (OjE) such that lim„_j.oo £n = there exists a subsequence {en^)ken cind a family of 
Borel probability measures on Zq , {pt)teTSi , satisfying 

M{ge„^{t), keN) = {fit}, 

for any t S K . 

Proof. We only sketch the proof and essentially indicate the points which differ from [5', Propo- 
sition 3.3]. Let us write 

G,it,O^Ty[g,{t)Wm ■ 
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By using Proposition IB. II and (1 + N) < 2(1 + N^) < 2(1 + S'^), we can prove like in [B^ that 



for some constant C > . We have 

\GS,V)-G,{s,0\ < \Tr[rg,{t)-~g,{s))W{0]\ + \G,{s,0-G,{s,v)\ . 
On the other hand, by making use of Lemma 13.81 we know 

4 



(27) 



\TrMt)-g,{s)]Wim < 



f Tr[^,(u;)^£^-i6,(«;,Cr^'=]du; 
< Go\t^s\\\{l + S,)'^^g,{l + S,nc^ 



in) 



X sup \\il + S,)-^/^ 
we[t,s] 

< Cl|^-,s|(l + |ek)^ 



when ^ e ff^(R'^) . Taking an approximation ^„ e _ff^(R), n e N, such that hm„_>.oo ||^— = Oj 
2i = 7?^(R") , and taking the hmit as n —j' 00 of the left-hand side with the help of ([ST)) , allows 
first to extend the previous inequality to any ^ e -Zi . 
Thus, we conclude that 



\GS, r,) - GM 0\<c{\t~ s\mz, + If +\r,- ^Hl, + ICIlJ , 



(28) 



holds for all (s, ^), (t, 77) e M x Zi , uniformly w.r.t. e E (0, e) . Remember also the uniform estimate 
|G,(s,OI<l. 

Now, we apply the same Ascoli type argument used in [6l Proposition 3.3] in order to prove the 
existence of a subsequence (ent )fc and a continuous function G{., .) : R x Zi — C such that Ge^ {t, ^) 
converges to G{t, ^) for any t & R and ^ € Zi. Further from (^7)) we can extend G(., .) to a contin- 
uous function on M x Zq . An "^/3" -argument shows that for any {t, ^) e M x Zq, hm„_i.oo G^^ {t, ^) 
exists and equals G(t, ^), so that G(t, .) is a norm continuous normalized function of positive 
type. Therefore, for any < G M, G{t, .) is a characteristic function of weak distribution (or projec- 
tive family of probability measures) jit on Zq . Finally we end the proof as in [51 Proposition 3.3]. □ 



3.2.3 An additional compactness argument 

Here, the compactness assumptions (^^[3]) is translated into some compactness property of the Wick 
symbol 61 . It allows to refer indirectly to Proposition 13.31 and to take the limit as e — in the 
term with j = 1 in ([22]) . With the rough estimates used in Proposition 13. 9| the terms in (f22|) 
corresponding to j > 1 with a factor e^~^ will vanish as e — > . The next Lemma applied with 
Qeis) in the integral term of (|22l) . will end the proof of Proposition 13.61 

Lemma 3.10. Let gi^ be a family of normal states on % satisfying for some constant C > the 
estimate 

Tr [(1 + Se)Qe{l + S^)] < C uniformly w.r.t e G (0, e) . 
Here is given by Definition \3.5\ . Assume that Ai{Qe) = {iAj then for any ^ G Zi , 

limTr[£<,W^(Ofoi(s,0'^^''] = / e^*° <«-->6i(s,0N d^i{z) . 

Proof. The polynomial fci(s,^) G T'l 2 + "7^2,1 splits into two similar terms, namely 

Bi(z) = (e® and ^2(2) = (^^^ T4(z ® 0) ■ 



14 



with their associated operators 

Bi = ((ei «) 11) e £(L2(M2<i)^ L^{m.'^)) and B2 = ^jT^ (H «) |0) G 'C(i^(K'^), • 

Let X G Cg^W with x{x) 1 if |a;| < 1, x(a;) = if |a;| > 2 and < x < 1. For m e N* wc put 
Xm(a;) = x(^) and define 

^l,m =Xm(|i?.|)Si(ll®Xm(|C..|))52 and Ba^™ - ^2 (11 ® Xm ( I ^r. I ) ) ^sX™ ( I ^r. | ) 

as bounded operators in C{L'^,{R'^'^), L'^{M.'^)) and £(L2(R'^), L2(R2rf)) respectively. We claim that 
both operators Bi^m and -62, m are compact. Actually, B2,m = ^1 m and 

Si,„, = i(n(»e--^^2)((e-^^|,»]l)(e--i^-2;^„(|i?^J)y(a;2)x„,(|i^,J)e'-i^^^)(^^^^^ 
Moreover, the linear norm continuous application 

maps Hilbert-Schmidt operators into Hilbert-Schmidt since 

||((e'^^CI®ll)e-"^''-(l®^)||£2(i2(R2.),i2(R.)) = \^\zo\\A\\cHZo) 

comes by computing the Schwartz kernel with ||^||£2(^2(^)\ = / \K{x,y)\'^ dv[x)dv[y) . Hence it 
maps compact operators into compact operators, because the space of compact operators, is 
the norm closure of £^ in £ . Therefore, by taking A — Xm(|-Di:|) V{x) Xm{\Dx\) which is compact 
by assumption (^[3]), we conclude that Bi^„i and i?2,m are compact. 
Now, we write for j — 1,2 

\T^[g,W{OBY'^''] - Ai(e^*'= < |Tr[e,VF(0(Sf ^^'^ - (29) 

+ |Tr[g,W^(OSf^^'=] - M(e^*'= ^^^^ B,,n{z))\ (30) 

+ |^(e^*^ ™(z)) - /i(e^^«'= <«">Bj(^))| , (31) 

with Bj^m G 'P^gi^o)- The right hand side (PH)) converges to owing to Proposition 13.41 Since 
s - hm™_^ 

00 Xmd-Dxl) — 111 polynomials i?j ^(z) converge to Bj[z) for any z € Zq , while the 

estimate 

|S,;™(^)|<c|ekl|(l-A)-VV|||z||„ 
holds true uniformly in m for some constant c. In addition Tr[p£N'^/^] < C implies that 

/ kl|„rf/i<C. 

JZo 

Therefore by dominated convergence theorem the right hand side (|3ip tends to when m — >■ 00. 
It still remains to prove the convergence of the r.h.s of ([2^ . Indeed, writing 

x{S, + l)-i(Sf ^^'^ - B]^r'=)(5e + 1)"'] 
and then by using Proposition 12.41 we obtain the estimate 

mg^wiOiBY''" - Bf^^'^)]!! < c\\{s, + i)-i(sf-'= - i?]^r)(^e + ir'Wcm ■ 

By functional calculus of strongly commuting self-adjoint operators we see that (S,. + 1)^^(VN + 
dr(l — A) + 1) is uniformly bounded with respect to e e (0, e). Now, applying Lemma [6.41 (with 
A = 1 — A), we conclude that 

mg^WiOiBY'-' - Bf^^'^)]!! < 11(1 - A.J-i/2(4. _ 4.„J(i - A.J-V2|| 

< leUo 11(1 - A)-'/^V\\ciZo) 11(1 - A.)-V2(i _ Xm{\D.M\ciZo) 

Since by functional calculus ||(1 — Aa;)~^/^(1 — Xrwd^'il))!! ^ the r.h.s of (^9]) goes to when 
TO — > cx) uniformly in £ S (0,e). Finally, a "(5/3-argument" with the established convergence of 
(IMl), dSni) and dni) yields the resuh. □ 
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3.3 Asymptotic a priori estimates 

In this section, a priori information on Wigner measures are derived from a priori estimates on the 
state . In particular, we shall prove the next result. 



Proposition 3.11. Let be the operator given by Definition \3. 51 and assume that the family of 
normal states {Qe)^^^Q ^) satisfies 

ya e N, 3C„ > 0,Ve e (0,e), Tr [(1 + S,)g,{l + S,){1 + N)"] < C„ , 

and A4{gs, e € (0, e)) = {//} . Then the measure fj, is carried by Z\ , its restriction to Z\ is a Borel 
probability measure on Z\ and 

\A%M%, df^iz) = I \z\%M\o Mz) < +0O. (32) 

Zq J Zi 

We will need the two next Lemmas. 

Lemma 3.12. Let b be a non negative (self-adjoint) operator on V Zq and assume that the family 
of normal states {gs)e£{o,e)! with Tr [p^N"] < Cq for all a G N, satisfies 

Tr [pefo^'"''] < C and e G (0, e)) = M ■ 

Then Zq 3 z ^ b{z) = {z'^p , bz®^) G [0, +oo] is a Borel function on Zq and J^^ b{z) d^{z) < C . 

Proof. When b £ V^piZo) has a compact kernel b we know after Proposition (see [H Corollary 
6.14] for a complete proof) that 



C > lim Tr [g.b^'''''] = f b{z) dfi{z) . 

J Zo 



We use the fact that b }jV/ick jg operator monotone, in the following sense: if the (possibly 
unbounded) non negative operators 5i, 62 in Zq satisfy 62 > &i > 0, then 6^"^'^ > b^^'^'' > in 
H. ^ 

By taking b e Z), for 6 G Vp^p{Z), as the supremum of 6„ with 6„ compact, we obtain firstly 

for all n e N 



C > liminf Tr 



J £— >o L 



Secondly, the monotone convergence yields 



b„{z) dfi{z) . 



C > sup / (Z^P , bnZ'^P) d^i{z) = / b{z) dniz) . 
neN J Za Jzo 

When b is unbounded, it can be approximated by 6„ = e ^{\/^ ^a)^ for n > 1 . Set bn[z) = 

Kz'^P) . The function b{z) = {z'^p , bz'^P) = sup„gN6„(z) is a Borel function on Zq as a 
supremum of a sequence of continuous functions. The uniform estimate 

Tr [fed < Tr [g^b^'^^] < C 

with the result for 6„ e C{\J^ Zq) gives bn{z) dii{z) < C , for all neN*. Again by monotone 
convergence, we get 



b{z) dfi{z) ~ sup / bn{z) dii{z) < C . 

Zo new* J Zo 



□ 
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Lemma 3.13. Let A he a non negative, self-adjoint with domain VIA), operator in Zq . Assume 
that the family {ge)ee{o,e} satisfies the uniform estimate Tr [^j^N"] < Ca , for all a G N, and 
A^(£»e, £ G (0,£)) = {fi} ■ Then the following implication hold: 



(Ve € (0,£), Tr [gedT{A)] < C) / {z , Az) d^jL{z) < C 

\Jzo J 

(V£ e (0,£), Tr {Q,AV{Af \ <C) ^ (^J {z, Azf d[i{z) < , 
(V£ e (0,£), Tr [^7,d^(A)2N] <C) ^ (^J^{z, Azf\z\\, dii{z) < . 

In all the three cases, the measure fj, is carried by the form domain Q{A) of A . 
Proof. The first implication is a direct application of Lemma 13.121 applied with 

b{z) ^{z,Az} , b^A , 6'^*"'= = dr{A) . 
For the second one we use 

dr(A)2 = ((z«52 ^ ^ A)z^^))'^""' + edT{A^) > {{z^^ , {A ® . 
We apply Lemma 13.121 with 

b{z) = {z , Az)'^ and b^A^A. 

For the last one we notice that N — dr(l) and dr(A) commute so that 

dr(A)2N > N ((z«2 ^ ^ A)z^^})'^'"' . 

With N = (I^^I^q)^*'^'^! the composition formula of Proposition IBT2] (extended to an unbounded A) 
says that N {{z»^ , {A (g) A)z^^))^"'''' = b^"''' with 

b,{z) = \z\%^,{z,Azy' + 2e{z,Az}\ 

Hence we get 

dTiA^N > {\z\%^{z , Az)Y''\ 

and we apply Lemma 13.121 with 

b{z) ^ \z\%^{z , Azf , 6= i(ll(KiA(KiA + A(Kill(Xi^ + A«)A«)ll). 

For the last statement it suffices to notice that the integrand is infinite in the Borel subset of Zq , 
Zo\ Q{A) ^{z£Zo,{z,Az)^ +00} . □ 



Proof of Proposition [3TTTt With 5"^ = dr(l - A) + N^, while dr(l - A) and N commute, we 
know 

(1 + 5,)(1 + N) >dr(l- A). 
Hence Lemma [3.13l savs that the measure fi is carried by Q{1 — A) — Zi with 

/ \z\l^ dpi{z)= I \zW dy.{z)<C. 

Let us check that /i is a Borel measure on (Zi, | \zx^ ■ The tightness property is given by the above 
inequality. According to [4] [54] [58] , it suffices to check that 

with = (u, (1 - A)m)2o , 
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is a positive type function which is continuous w.r.t ^ restricted to any finite dimensional subspace 
of Zi . 

Consider the regularized version 

r -2i7rRc (^VC,z)2n f -2i7rRc (^t-5-z)z„ 



with A= (1 — A) . For all ^ G Zi the pointwise convergence 

Vz e Zi , lim (— — 4^ , 2)^0 = , z)^^ 

n 

and the uniform bound 

|e I < 1 

imply the pointwise convergence of the integrals 

V^eZi, hm Gi.„(0=Gi(0. 

n— f oo 

But Gi.n{C) equals G((l + where G is the characteristic function of ji in Zq: 



2o 



Hence for every n e N, the function Gi_„(^) is a positive type function. As a pointwise limit of 
Gi.n, the function Gi is also a positive type function. 
For the continuity, the equality 



implies 

|Gi(e) - Gi(e')l < Mi - riz, ^ \z\z, di,{z) < 27r 1 + |z||^ d^x{z)^ le - C'lz, , 

and the function Gi is a Lipschitz function on Zi . This finishes the proof that is a Borel 
probability measure on Zi . 

For the inequality (I32p . it suffices to notice the inequality of (commuting) operators 

(1 + S,f{l + N) > (dr(l - A)2)N . 

Applying Lemma [3.131 vields 

/ |z||Jz||^ dMz)<G. 

□ 



3.4 Uniqueness of the mean field dynamics via measure transportation 
technique 

Now we are in position to prove Theorem 11.11 This will be done in three steps: 1) Writing a 
transport equation, in a weak sense in Zi for /i^ ; 2) Solving this equation as = ^(i, 0)*/^o when 
the initial state fulfills strong decay estimates; 3) Relaxing the strong decay estimates. 
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3.5 The transport equation on 

We shall need similar notions about cylindrical functions, as the one used in Zq and recalled in 
Appendix IB. 11 Let Pi denote the set of all finite rank orthogonal projections on Zi and for a given 
p e Pi let Lp^i{dz) denote the Lebesgue measure on the finite dimensional subspace pZi, with 
volume 1 for a Zi-orthonormal hypercube. A function / : Zi — > C is said cylindrical if there exists 
p € Pi and a function g on pZi such that f{z) ~ g(pz), for all z € Zi . In this case we say that 
/ is based on the subspace pZi. The set of (resp. S) cylindrical functions on Zi, is denoted 
by C^^yi{Zi) (resp. Scyi{Zi)). We shall also need C^^yi{Zi x M) , in which the algebraic tensor 

alg 

product C^^yi{Zi) <Si C^(M.) is dense. Finally the Fourier transform of elements of Scyi{Zi) is 
given by 

f{z)= [ J-i[/](C) e2-*^<= WAdO- 

According to the notations of [3_ and Appendix [Cl we consider the space Prob2(Zi) of Borel 
probability measures fi such that 

/ \z\l^ dfi{z) < +00 . 

JZi 

On this space, we introduce the Wasserstcin distance 

W^{pi,t^2)= inf / \zi ~ Z2\%^ d^i{zi, Z2) (33) 
p.er{p,u^i2) Jzl 

where V{^i,^2) is the set of Borel probability measures /i on Zi x Zi such that the marginals 
(IIi)*^ = ^.i and (112)*^ = ^2 ■ 

Proposition 3.14. Let and Qe{t) he the operators given by Definition \3.5\ and (j20p . Assume 
that the family of normal states (£'e)^g(o e) satisfies 

ya e N, 3C„ > 0, Ve e (0,£), Tr [(1 + 5^)^(1 + S,)il + N)"] < C„ , 

and consider a subsequence {sk)k<£N, £fe '^^^ such that 

M{ge,{t), fcGN) 

according to Proposition \ 3.6\ Then the measure fit is a Borel probability measure on Z\ which 
satisfies 

• J^^ \z\z-^\z\%g dflt{z) < C for some C independent 0/ t G M . 

• It is a Lipschitz function oft e R m Prob2{Zi) endowed with the Wasserstein distance (j33|) . 

• It is a weak solution to the Liouville equation 

dtfi + i {Vs,fit} = , 
in the sense that for f € 0(^1 ^ ^) 

f f {dtf + t{VsJ})dfit{x)dt^O (34) 
Jm Jzi 

Proof. The Proposition 12.31 as well as the commutations [e**^^,N] = [e**^e,N] ensure 
Va e N, 3C'^ >0,yte R, Ve G (0,£) , Tr [(1 + 5e)&(0(l + Se)il + N)"] < C'^ . 
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The ProDOsition l3 . 1 1 1 for any i e M, provides the first results, which also contain fl{t) £ Prob2(Zi) . 
It remains to check the continuity with respect to the Wasserstein distance W'^ and the Liouville 
equation. Replace ^ e Zq, by V^n^ in the equation (PT|) in Proposition 13.61 and integrate with 
T[m) dL^{z) with 

/(z)= / J-[/](z)e^-^^<^^«> dL^{z) 
when / g Scyi[ZQ) based on pZo . This gives for any t, e M and all / G 5cy/(Zo): 

f{z) dfit'(z)^ f fiz) dfit(z)+t f I {VsJ}{z)dfLs{z)ds. 



' Zq J t J Zq 

With the uniform estimate 

|6 j~ ^ \ U|4 



dUz)< / |z||Jz|i„ rf/I.(z)<Co, 

Zq J Zq 

while Vs is quartic, it can be extended by dominated convergence to any cylindrical polynomial / 
of total degree less than 4 . 

Consider now an increasing sequence of finite rank positive operators An on Zq such that 

supAi = (1- A). 

riGN 

Since jit and flf belong to Prob2 {Zi ) and are Borel measures on Zq carried by Zi , we have 

{jit' , fit) < / \zi - Z2\%^ dp.t{zi)djlt,{z2) = / \zi- Z2\%^ djlt{zi)djlt'{z2) 

and the monotone convergence ensures 



W\ft',ilt) < sup [ ( 



Zi- Z2, An{zi - Z2)) djlt{zi)djlt'{z2) ■ 



For a given n e N, and a fixed zi e Zq, the function Z2 — >■ fnizi, Z2) = {zi — Z2 , An{zi — Z2)) is a 
cylindrical polynomial of total degree 2 . Hence we get 



{zi - Z2 , Anizi - Z2)) dflt{zi)dilt'iz2) < / ds \{Vs, fn{zi, .)}\ djls{z2)djlt{zi) . 

Za-KZa •'[t,t'\ JZoxZo 

We compute 

{V, , /„(zi, .)} (Z2) = 4z Im (zf 2 , v;(;22 ^ (A„Z2))) - 4i Im {zf^ , V;(z2 (A„zi))) . 

1 /2 

From P7)) and with uj^o < l^lzi, we deduce 

|{V; , /„(Z1, .)} (Z2)| < C(|Ay2^2|zo + \AIJ^Z,\ZJ\Z2\%\Z2\Z0 <C{\Z2\Z, + kllzj k2||,k2|2o • 

The uniform estimate J^^ \z\'z-^\zW^ djir{z) < Co valid for all times, thus leads to 

w\^it.,^lt)<c'\t' -t\. 

In order to prove the last point, let us come back to the equation (PTjl . Again the estimate ([T7|) 
implies 

|im (z«, v;.(c®z))| < q^llJ^lzoieU-MR")- 

Since the measure fit is carried by Zi for all times with moment estimates, the first line of (|2ip 
shows that it still makes sense for ^ G . By taking ^ = ^/2tt{1 — A)?7 with 77 e Zi, we 

obtain 
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By integrating with J-i{g)(r]) dLp^i{z) for a cylindrical function g e Scyi{Zi) based on pZi, we 
get ^ 

/ g{z) djitiz) ^ I g{z) dfloiz) + i / , 5} (z) dfls{z) ds . 

J Zi J Zi JQ J Zi 

Hence for any / e Scyi{Zi), the function Ig : t ^ J^^ g{z) dp,t{z) belongs to C^(M) with 

dtlg{t)=i( {Vt,g)[z)dlit{z)- 



By multiplying the above relation by f{t)^ with ip e C(5"(R), and integrating by part proves ([34|) 
when /(t, z) = ip{t)g{z) . We conclude by the density oiC^^yi{Zi) Cg°(K) in C^^yi{Zi x M) . □ 



3.6 Uniqueness of the measure for regular initial data 

Proposition 3.15. Let and Qi^it) be the operators given by Definition \3.5\ and (|20p. Assume 
that the family of normal states (ge)jg(o e) satisfies 

VaeN,3a >0,Vee (0,e), Tr [(1 + + + N)"] < 

and M{ge,e e {0,e)) — {fio} ■ 

Then for any time t E M, the family {gs{t) = e^^'^' (>e'^^^^^)Ee{o,e) admits a unique Wigner 
measure fit = 0)*/io, where $ is the Hartree flow defined by ^ on Zi . It is a Borel proba- 
bility measure on Z\ which belongs to Prob2{Zi), which is Lipschitz continuous w.r.t time for the 
Wasserstein distance W'^ and which satisfies 



VteK, / \z\%^\z\l^dfit{z)<C. 
Jz^ 



Proof. We still start with the state Qeit) defined in (|20p . Proposition 12.61 says that the group $ 
associated with ^ and the dynamical system $ associated with 

idtz^v{t,z) , v{t,z) = e-''^{[V '^\e'"^z\^]e'''^z) 

are well defined on Zi . Further it gives the estimate for the velocity field 

<|ini-A)-i/^|||z||Jz|zo. 

When fit is the Wigner measure defined for all times and associated with a subsequence {en^)keN, 
we obtain 

VteR, / \v{t,z)\l^dfit{z)<C [ \z\%^\z\%^ dfitiz) < C . 

JZi JZi 

With Proposition 13.141 this means that (/if) fulfills all the conditions of Proposition IC.8I and we 
deduce that fit — 0)*/xo . But this uniqueness implies A4{ge{t),e G (0, e)) = {fit} for the whole 
family {ge{t))£e{o,e) ^nd all times t e M. 

Going back to g^{t) — e~^i^" gs{t)e^i^' , it gives /xt = $(t,0)*/.to . 

The additional properties are the ones coming from Proposition 13.141 □ 



3.7 Evolution of the Wigner measure for general data 

We follow the truncation scheme used in j6j. When the initial data satisfies only 

\\{N + H°y/^g,{N + H^)^/^ <Cs 
for some S > 0, we approximate by 
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with X-rI": ^) = xItI I tI) with i? +00 and < x < Ij X G C5"(M^) and x = 1 in a neighborhood 
of . The time evolved state is defined by 

The assumptions ensure that for ah times 

with ly independent of {t,e) and hm/^^oo i^iR) = 0. We recall the Proposition 2.10 of [6]. 

Proposition 3.16. Let i0l)ee{o,e)' ~ families (or sequences) of normal states on % 

such that Tr [gi^ N'^] < Cs uniformly w.r.t. e £ (0,£) for some S > and Cs & (0, +00). Assume 
further M{gi,e G (0,£)) = {^j} for j = 1,2. Then 



J Imi - Ai2| < lini^inf Wei - glWcHn) 



Proof of Proposition 3.16, end of the proof of Theorem ll.lt For R G (0,+cx)), the state 
Qe^r fulfills the conditions of Proposition l3 . 1 Sl except the uniqueness of the Wigner measure at time 
t = . Out of any sequence (en)nGN, a subsequence (erifc)fceN can be extracted in order to ensure 

Thus after this extraction we obtain 

VteM, X(&,.,.K(t),fceN) -{$(t,0),Mo,j?,} • 
Take t e M and let fi belong to M{geit),e G (0,e)) . There exists a sequence (e„)„gN such that 

A^(e,„(i),neN) = {/^} . 
After extracting a subsequence like above and by using Proposition 13. 16"! we obtain 

J |/i-$(i,0)*Mo| < J |Ai-$(i,0)*Mo,fl| + J \mM~ ^^o\<MR)■ 
Taking the limit as i? — > cx) implies /i = <i>(t, 0)*/io and therefore 

X(&,£e (0,e)) = {$(t,0)*Aio} • 

This also proves that limi^^oo Jz l/^t ~ ^(t^ (^)*IJ'0,r\ = , while all the measures 0)*/io,_R are 
Borel probability measures carried by, and on, Zi . This implies that is carried by Zi and is 
also a Borel measure on Zi . This ends the proof of Theorem 11.11 □ 



4 Complements 

Additional results are given in the three first paragraphs, concerned with the BBGKY hierarchy 
or the propagation of energy. The fourth one shows some examples and the last one is an informal 
discussion about the classical mean field problem. 

4.1 BBGKY hierarchy 

Although the analysis here is different from our previous work 6 it is possible to combine them, 
in order to strengthen the result of Theorem ll.il It is also interesting to reformulate our result in 
terms of reduced density matrices since, in the literature, several mathematical results on mean 
field limit use the BBGKY hierarchy method ([8l[9l|42] for example). For a family of normal states 
{9e)ee{Q.e) ^ ^lud r G N, the reduced density matrices 7!*'^ € £^(L^(E'*p)) is defined according 
to 



Tr 

with the convention that the right-hand side is when Tr [pedzp'')^*'^'^] — 
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Proposition 4.1. Let {ge)eeio.e) family of normal states on %, satisfying the hypothesis of 
Theorem ] 1.1[ with a single Wigner measure fiQ such that 

Va e N, lim Tr[6ieN"] = / jzp" d^oiz) < +00 . (36) 

Jzo 

Then for all t G K, the convergence 

limTr[e,(t)6^'"'=] = / 0)z) d/io(^) = / b{z) dfit{z) 



holds for any b G Vaig{Zo) = ^l^len'^p^yi^o) , with fit = $(t,0)*/io • 
Finally, the convergence of the reduced density matrices 

lim 7i^)(i) = I [ |z«0(z«^'| df,,iz) , 

Jzo 1^1^ '^l^tiz) Jzo 

holds m the C^{Ll{W^P))-norm for all p e N . 

Proof. We know from Theorem 11.11 that the the family of normal states (gie(i))£g(o,£) admits 
a single Wigner measure fit equal to $(i,0)^/io • Since the quantum and classical flows preserve 
the total number, the state ge{t) satisfies as well the condition p6p for any time i € R. The [6j 
Proposition 2.11, 2.13] provides the claimed results. □ 



4.2 Moment upper bounds 

In f4^, it was proved that the sole a priori estimate Tr [psN'^] < Cs for a given S > (possibly 
small), with M{g£,e € (0,e)) = {fi} leads to 



Izp'^ dfi{z) < +00 . 



\2S 

1^ 

The a priori estimate, assumed in Theorem 1 1.1 1 at time t — 0, leads to 
according to the following result which is a variation of Lemma 13.121 

Proposition 4.2. Let {A,T>{A)) be a self-adjoint operator on Zq such that A > I . If the family 
of normal states {ge)e(£(o,6) satisfies Tr [gie(dr(A))*] < Cs for some 6 > and M. {qs, £ G (0, e)} = 
{/i}, then 

{z , Az}^ dfi{z) < +00 . 
Proof. By Wick calculus (see Proposition IB. 21 when A is bounded), one gets 

dr(.l)'= > ((z, Az)'')'^'"' ,VfceN. 

Let (ej)jgN be an orthonormal basis of Zq such that ej £ 'T>{A) for all j £ N, and set 

Aj = X:^'/^|e.)(e.|Ai/^ ^ y: \A^/'e,){A^/'e,\ . 
j=o j=o 

The inequality A'^'^ > A'^^ holds for all J S N, while b — >■ fe^'^'^ is operator monotone when 
restricted to operators b acting in V^*^ Zo . Therefore, we obtain 

/ „ \ Wick 

\k \ 



(i+dr(A))" = j2^ndm)'^>[J2c^^{z,Azy 

> lj2c'n{z,Ajz)A ^[{l + {z,Ajz)rr'^\ (37) 
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We now use the same argument as the one used in [3] when A = Id, relying on the semiclassical 
calculus in finite dimension (see [12 ) |40j [ 50] | 52 ) |56| ) . 

Let pj be the orthogonal projection from Zq onto ©j^o^^^^^i • '^^^ symbol 1 + (z , Ajz) is a 
cylindrical symbol based on pjZo ■ Since ker A^/^ — {0} and eo, . . . , ej are linearly independent, 
the symbol 



{l + {z,Ajz)) = l + J2\{A' 

3=0 



is an elliptic symbol on PjZq 



-<j+i 



in the Hormander class 5(1 + |z|?,j+i. 



\dz\ 



The 



(38) 



(39) 



functional calculus of Weyl e-quantized elliptic operators in finite dimensions gives 
Vs e R, [(1 + {z , Ajz))'^^y'Y > (1 - Cj,,e) [(1 + {z , Ajz)y]'^'y' . 

The finite dimensional comparison of Wick and Weyl quantization, also gives 

Vn e N, [(1 + {z , Ajz))^]"^''' > (1 - Q„£) [(1 + {z , Ajz))'^^^']" . 

From ([37l)(p9l) and the operator monotonicity of _B — ^ _B* for t e (0, 1] , we deduce 
Vs e R , (1 + dr{A)r > (1 - C:,,,e) [(1 + {z , A,;z))^^^'] ' , 

and ([38]) gives 

vs e M , (1 + dr{A)y > (1 - c:;,,£) [(1 + {z , A^z))^]^^^' . 

The definition of Wigner measures, recalled in Theorem 13.11 says 

limTr = /" b{z)dfiiz), 

J Zq 

for all b £ ScyiiZo), in particular the &'s based on PjZq . Take now s = (5 in (|40p . The a priori 
estimate 



(40) 



Tr 



f?, [(l + (z, A/^))']'^'''' <{l + Cj^se)T^[Q,{l + dT{A)Y] < Cs{l + Cj,se) . 



and the ellipticity of (1 + (z, Ajz)Y allows to extend the above convergence to any cylindrical 

|2 \<5 l^^fl^J+i 
•IcJ+i 



6 = / o pj with / e S{{1 + |z|?.j+i)'', f ?2 ) . In particular, this leads to 



C^- > limsupTr [^,(1 + dr(v4))*] > lim Tr [(1 + (z , A.jz)f]^"'^] = I {l + {z,Ajz)Y dfi{z) 



Since A — supj Aj with Aj' > Aj for J' > J, the monotone convergence implies 
Cs> [ {1 + {z , Az)^ dfi{z) > [ {z , AzY d^i{z) . 



□ 



Proposition 4.3. Within the framework of Theorem \l.l\ with the assumption Tr[(N + _ff^)'^gie] < 
Cs for 5 < Q , the measure satisfies the additional estimate 




\z\l dfitiz) < Cs 



for all times t G R . 
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Proof. The functional calculus of commuting operators implies 

csil + < (1 + N + H'^J < (1 + SeY 



Thus the initial state , satisfies 



Tr 



&(1 + SeY'^ 



From Proposition [231 we deduce 

Since i? — > i?'' is operator monotone for s e (0, 1], this implies 

e'i"' (1 + S.y/^e-'i"^ < + S,)'/^ 

as soon as 4 < 2 . The inequality 



Tr 



fe(0(dr(l - A))*V3l < Tr L(i)(l + 



and the previous Proposition 14.21 applied with A — (1 — A) , yields the result. □ 
A more accurate version of this last result is given below by making use of the conservation of 
energy. 

4.3 Convergence of moments and energy conservation 

For a family (pe)gg(o of normal states with a single Wigner measure /ig the condition p6p is an 
important and non trivial assumption. Indeed, we proved in W the following equivalence 

Va e N, Imi Tr [g,N"] = £ |z|2" dM^)^ ^ (^V6 G Vaig{Zo), lim Tr = ^ b{z) dfio 



Hence the condition ([36]), although it involves only the number operator, is exactly the one which 
leads to a good asymptotic behaviour of the reduced density matrices. 

Proposition 4.4. Let ige)ee(o.6) family of normal states on %, satisfying the hypothesis of 
Theorem ] 1. 11 with a single Wigner measure fiQ. Assume Tr[£i£ N"] < Cq, uniformly in e € (0,e) , 
for all a eN . 

Then for every a G N, the quantity 



liminf Tr[^.,(^)N«]- / |z||^ dfit{z) 



does not depend on time when Qe{t) — e^^i^^ Qec'^c^^ and A4{ge{t), e G (0,£)) = {/it} . 

The condition p6p is satisfied by (f?£(i))Ee(o,£) ^t^'^ f^t, for all times f G as soon as it is true 

for one to £ R . 

Proof. According to Theorem ll.il we know that A4{ge{t),e G (0, e)) — {fit} with fit — <&(i, O)*fio ■ 
The conservation of the [.j^Q-norm by the nonlinear flow 0) yields 

Jzo 

for any i G M. On the other hand, and N are strongly commuting self-adjoint operators and 
therefore Tr[fe(t) N"] = Tr[£ieN"] for every a G N . □ 
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Proposition 4.5. Let {ge)eeio.e) a family of normal states onH with a single Wigner measure 
fiQ satisfying the hypothesis of Theorem \1.1\ with S — 2 and condition (1361) . Then for any i € R 



lim Trfpgft) iJe] = / h{z,z) d^i 
Jzi 



(z) e (-00,00) (42) 



where Qe{t) = e ^^i^' g^e^*'^' , fit — 'I'(^, 0)*/xo o,iT-d h{z^z) is the classical energy given in and 
both sides of the identity do not depend on time. 

Proof. With the energy conservation, it sufBccs to prove ([^^ for i = 0. Let x ^ C'o^C^) such that 
< X < 1, x(s) = 1 if |s| < 1 and x{s) - if |s| > 2. For m e N* , set Xm{x) = x(^) . Let Bi{z) 
and B2{z) be respectively the polynomial (z, —Az) and B2{z) — V{z) = ^(z*^^, V{x — y)z^'^) well 
defined for z G Zi. Remember that although the kernels of Bi and B2 are unbounded operators 
their Wick quantization still have a meaning as densely defined operators on H (see Appendix |B|) . 
We write for j = 1 , 2 

\Tr[g,Bp^'']- [ B,iz)dM^)\ < |Tr[e, (^f - -'=)] | (43) 

JZi 

+ |Tr[e,i3]^-'=]- / S,,,„(z)d/io(z)| (44) 
Bj,m{z)dfj.o{z) - Bj{z)dfio{z)\, (45) 

JZi 

where Bi,™(z) = (z, -Axm(-A)z) and B2,m{z) ^ ^{z'^'^ ,V{x - y)x,n{~A.^)z'^^} . We observe 
that Lemma FB. 61 leads to 



||(dr(-A)+N+i)-Mr(A(i-x™)(-A))(dr(-A)+N+i)-i|| < ||^-_^(i-x™)(-A)|| ™4°°o 

and 

||(dr(-A) + n2 + i)-i/2 ^B^^ck _ QWu^k^ (dr(-A) + + i)-i/2|| 

< Cv\\{l - A)-i/2(l - x™)(-A)j| '"4°° . 

Therefore the r.h.s (|43p tends to when m ^ 00 thanks to the regularity of g^ and by noticing 
that (dr(-A) + N2 + l)i/2(dr(-A) + N + 1)-^ is bounded. Now, since Bj,™, J = 1, 2 belong to 
VaigiZo) then by the statement (PT|) . proved in [BJ Proposition 2.12], the r.h.s ([^^ converges to 
when e — >■ 0. Further, by the dominated convergence theorem and with the help of Lemma 13.131 
we see that the r.h.s (I45p vanishes too as m — >■ 00. Hence an (5/3-argument gives 

lini Tr [g, Bf'"'"] = j Bj (z) dfia (z) for j = 1 , 2 . 
Thus we have proved (|42|). □ 

4.4 Examples 

We give here two examples, other can be seen in our previous articles [4] [5] [6]. The first one recalls 
that the transport of the Wigner measure takes into account some correlations. The second one is 
about the mean field dynamics of states, which do not satisfy (I36p and makes a connection with 
Bose-Einstein condensation. 

4.5 Deformed torii 

For two elements ipi,ip2 € Zi C Zq such that — \\ip2\\ ~ 1 and {ipi , ?/'2) — 0, the space Zq 

can be decomposed into 
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This decomposition is second-quantized into the Hilbert tensor product 

n = r,(Zo) - r,(c^i) ® t,{c^2) ® r,(V'^) , 

which ahows an analysis by separating the variables. The number observable is now 
N = (Ni ® 11 (g) 11) ® (11 N2 (g) 11) ® (11 (g) 11 ® N') , 

simply written as N = Ni + N2 + N' and where Ni, N2 and N' are respectively the number 
operators on rs(C'0i), rs(C-02) and Fs (■(/'''") ■ Consider in this decomposition, the state 

Q,^Bl®gl®{\n'){n'\) 

where \^') is the vacuum state of Ts{ip'^) and 

with lim erii = lim en2 = — ■ 

e-i-O £^0 2 

In H = TsiZo), this state is explicitly written (see |6]) as 

= |^V(„i,n2)^^^V(ni,n.)| (4g) 

ni times 712 times 

with .,1 , , a*(^i).^a*(V^i)V(^2)...a*(^2)|»). (47) 

The state satisfies 

limTr[NVe] =(^ + ^)'= = l, 

owing to N = Ni + N2 + N' . Moreover, with gB]) gZl) , N + = dr(l - A) and the help of Wick 
calculus, it also fulfills 

lim Tr [(N + H^)g,\ = . 

Meanwhile the separation of variables allows to compute explicitly the (it is unique) Wigner mea- 
sure of {ge)ee(0,6) 

/-io = (5^, <8<5^, «)<5o on Zi = (CV^i) X (C?/'2) X , 

1 1 /■^'^ 
with (5f ' =^J^ ^e^'u d0 . 

We get 

/ Izp'^ dfioiz) = / (kiP + \Z2\' + k?)' dfioiz) = 1 = lim Tr [N'^g,] . 
Jzi Jzi 

Hence all the assumptions of Theorem 11.11 and Proposition 14.11 are fulfilled. 

This measure is carried by a torus in Zi better described by using an other orthonormal basis of 
CV-i ® CV'2: 

\/2 ^/2 

V'o = — (-01 + V'2) , V'f = i — (V'l - ■02) , 

0^ = cos{(p)'ipa + shi{ip)ip^ , 
^(e'Vi+e^«'02) = e^'^V'^, 



1 



2ir 
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Two elements e^^^jj^p and e*^ V^s' '^^ the support of /io are equal when 

(6*' = 6 and = ip) or (6*' = 6* + tt and i^' = + tt) . 

Hence a one to one parametrization of the torus can be done hy £ [0, 27r) and (z [ip,ip + n) . 
Let i^ipit) — ^{t, 0)ip^, be the solution to the Hartree equation 

\ ^p^{t = 0) = ip^ ^ cos(93)-0i + e~*"i sin((/?)'(/'2 ' 

The gauge invariance of the equation says that for any 9 S [0, 27r], e^^ip^{t) = 0) [e^^Vv'] ■ By 
applying the result of Theorem 1 1.11 and Proposition 14. 1 1 we get 

^*^2^Jo '^^^A^ la Jo '^^'^^ 

VpeN, lim7i^)W = ^£^|[V^^(irO([V'^(tr''| dp. 

Since the Hartree flow is non linear, the complete hierarchy of reduced density matrices have to 
be taken into account if one wants to write evolution equation for them. More simply, they can 
be computed after solving an autonomous equation for the Wigner measure. Due to the nonlinear 
term the dynamics of correlations is by far nontrivial. This can also be thought geometrically: 
The initial measure is initially supported by a torus which lies in a 2-dimensional complex vector 
space (think of the circle in the plane R^'offiKV'f ); along the time evolution, the measure fit is still 
carried by a torus in Zi , which nevertheless, is a priori not embedded in any finite dimensional 
subspace . 




Fig.l: Evolution of the measure initially carried by a torus in C?Ao © . 
The complex gauge parameter e*^ is represented by the small circle. 

In Figure 1, the deformed torus for time t ^ 0, has to be imagined in the infinite-dimensional 
phase-space Zi G Zq ■ Contrary to the picture, there might be no intersection with the real plane 

This discussion can also be extended to higher dimensional torii after taking a finite (or countable) 
orthornormal family {ipn)i<n<N for building the initial states Qe with a measure Y[j=i ^\-^- (^^^ 
i)- 

4.6 Propagation without the convergence of moments 

In [3] we considered the thermodynamic limit of a free Bose gas on a torus with the one particle 
energy given by — A . We showed that when there is Bose condensation the condition ([55)1 fails 
and illustrates what we called a dimensional defect of compactness, in opposition to the phase 
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space or microlocal defect of compactness (see [1H][SS])- Others examples were given. In [5] the 
propagation result for bounded interactions but without any compactness condition, cannot be 
applied for such initial states. With Theorem 11.11 the propagation holds for these kind of initial 
states. Since our analysis is valid on R'' the analysis for the torus does not apply directly and we 
adapt the presentation of the Bose-Einstein condensation. 

Moreover the dimensional defect of compactness which plays with all the directions of the phase- 
space Zq — L^(R^), can be geometrically thought in the one particle phase-space T*]R'^ . The 
condition which leads to estimates of d^, suggests that the dimensional defect of 

compactness is due to mass going to oo in the position variable rather than in the momentum vari- 
able, in T*W'^ . The mean field limit that we consider here, can be tested by using the harmonic 
oscillator Hamiltonian A — —d^ + ^ — j . The motivated reader will then see that the dimensional 
defect of compactness is incompatible with the condition 

We work in dimension d> 2 . Let eo be an L^-normalized C(^{M.'^) supported in the hypercube 
(—5, ^)'^ and set 

VfceN'', efe(a;) = eo(x-fc). 

The family (efc)fegN is orthonormal in Zq — _L^(R'') . The spanned Hilbert subspace and the 
corresponding orthogonal projection are respectively denoted by Ze and He , Ran He — Ze . Note 
that 

ne(-A)ne= (^l^eop) He. 

Consider now the self-adjoint operator defined on Zq = L'^{R'^) by 

d 

^ = ^ \k\\ek}{ek\, |fc| =^fcj, 
kem j=i 

which restricted to Z^ is unitarily equivalent to the harmonic oscillator Hamiltonian A = —d^ + 
^ — I on R . We use the tensor decomposition 

Zo^Ze®z^ , n^r,iZo)^r,iZe)®T,{z^), 

N = Ne® 11+]1®N^ =NeH-N^ , \n) ^\ne}<»\ne) , 
VBe e CiZe) , \\Be\\ < 1, r{Be) ^ mi){^i\) = TiUeBeHe) , 

VB e c{Zo) , ||s|| < 1, TiUeBne) - r(ne o b\zJ <g> {\n^){n^\) . 

In particular the last relation with B = e*^*"^ differentiated at time t = gives 

r(ne)dr(-A)r(ne) = (^^^ |Veop^ . 

'tensorized) quasi-fre 

r(nee-^^('^-'^-)ne) (48) 



Consider on H, the e-depcndent gauge invariant (tensorized) quasi- free state 

1 



Tr [r(nee-''-(^-^«)ne)] 



^ T{UeZ,e-f'^^Ue) ., 



-TrlTilieZ.e-P^^Iie)] 



^ -r(ne)r(Zee"^'^i^<='^^(-*' )r(ne) . 



The chemical potential /^^ is negative of order e^~^^'^ and the temperature is large according to 



Z^ = = 1 - — , P,^ e^/'^ . 

vc 

With the e-dependent definition of a(/), a*(/), [0(5), a*(/)] = e{g,f), and W{f), this quasi-free 
state is characterized by the two-point function 

Tv[g,a*if)aig))] ^ s {Ueg .Z.e-f'^^^il - Z.e-f'^^^r'Uef) . (49) 
or Tr[e,iy(/)] =cxp[-£(ne/, (l + Z,e-^--4=)(l-Z,e-^=^-)-ine/)/4] . (50) 
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In particular the total number (multiplied by e) is given by 
Tr[geN] = Tr [f^.N^] ^ Tr [g.a* (efe)a(efc)] 

= vc + >^ + T'{e) with lim r{e) — , 



and 



du^\S'^-^\ 



1 - e-l"l 



+00 g-t 



1 - e- 



f^-^ dt, {d>2). 



(51) 
(52) 



We deduce 



lim Tr [p^N] = vc 



hm Tr = lim Tr [e,r(ne)dr(-A)r(ne)] - f / iVeop) 



and the condition (|4]) of Theorem II .11 is satisfied. 

Actually i^c > corresponds, in the analysis of the free Bose gas (see [4]), to the density associated 
with the condensate phase. In the scaling that we consider, it is the other part which produces 
the dimensional defect of compactness. Let us compute the Wigner measure, by considering the 
limit of Tr [geW{V2Trf)] as e . With 



/ = E + ' I/I' = E IZ-^I' + l/^l' - Ine/P + 1/ 



±\2 



the expression (f50|) gives 



Tr 



g''W{V27rf) 



= e--'|n=/lV2 X exp 



—en 



Ei/^ 



{1- Z.e-^'^-'W) 



The family {ge)ee(o.e) admits the unique Wigner measure 



(53) 



Mo 



> So on Zq — (Ceo) x 



which is carried by Ceo C Zi and which can also be written 

"+030 



Mo 



e 



(dz) 



(SV du. 



/Ceo T^i^C 

In particular, we get 

/ \z\%^^dfio{z) = i^c < i^c + = lim Tr [p^N] 

and ^ |z||^rfAio(2) = ^^c (^^jVeop^ < (i^c + ^^) (^^jVeop^ = Ihn Tr [e,dr(l - A)] , 

and the condition ([55)) does not hold. Even at time t = 0, no formula is available for the reduced 
density matrices in terms of the Wigner measure. Nevertheless the time-dependent Wigner measure 
of Qsit) = e 



~i-He 



g^e^e^'^ is given by Theorem II. 1[ since the condition 1^ is verified. Consider the 



solutions to the Hartree initial value problems 

idt^pu = -Aipu + {V * \ipu\'^)^u 
ipuit = 0) = ^eo , u € (0, +00) . 

Then the Wigner measure of gs{t) = e~^i^^ g^e^i^^ is given by 



Mt 



+ C30 — — — 

V(t) du. 



Again like in the example of the previous section, the measure /it is carried by surface containing 
and topologically equivalent to C, but this 2-dimensional surface does not remain a priori in any 
finite dimensional subspace of Zi for t ^ Q . 
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4.7 About the classical mean field problem 

The classical analogue of our analysis is the derivation of the Vlasov equation 

dtf + v.d^f - ^{d^,Vfix,t)).d,f = 
f{t,x,v) = fo{x,v) 

Vf{x,t) * gf{x,t) , gf{x,t) ^ J^^ f{x,v,t) dv 

where f{x, v, t) represents the particle density in the 1-particle phase space R^"*^ , from the classical 
Hamilton many body system 

I ; 

\ i^--,iw(J:f=id..V{x,-x,)) ' * = 

in the limit N ^ oo . This problems is still open for singular potential and C. Villani, in a recent 
survey article about the Landau damping 62: quotes the work of Hauray-Jabin j36| as the most 
advanced one in this direction. It works for a potential such that \WV\ ~ 0{\x\~''), s e (0, 1), and 
does not include the Coulomb interaction. 

Indirectly our result, justifies the mean field model up to Coulomb interaction in dimension d = 3 . 
In and more recently [2] , the Vlasov equation is proved to be the semiclassical limit of the semi- 
classical Hartree equation. This means that there are two "semiclassical" limits, one in the phase- 
space L^(R'';C) with the small parameter another one on the phase-space T*R'^ - R^'* for 
the one particle nonlinear problem. This double asymptotic regime is well presented in [25[ \27\ 135) . 

A possible strategy, for deriving directly the classical mean field limit from the classical many 
body problem, consists in adapting our approach by, as usual, replacing traces by integrals. For 
information, we refer the reader to the presentation [18] by J. Derezinski of the classical analogue 
of second quantization. Of course classical mechanics, although living in the commutative world, is 
often more singular than quantum mechanics, from the analysis point of view. With the Coulomb 
interaction, the Kustaanheimo-Stiefel desingularization of the hamiltonian fiow may be useful (see 
a.e. [TSllMlllTJIlSllig). 



Appendix 

A Commuting self-adjoint operators on a graded Hilbert 
space 

We briefiy study the general structure of self-adjoint operators on a graded Hilbert space. Prop- 
erties collected in this section are useful for the analysis of the quantum Hamiltonian ([T]). In this 
appendix, the small parameter is not required and we work with e = 1 . 

Remember that a graded Hilbert space "H is a direct sum of Hilbert spaces H„, n e N, of the form 

oo 

Let (A„)„gN be a sequence of self-adjoint operators where each A„ acts on Tin- We define the 
operator 

oo oo 

2?(A) = e -H : ^ ||A„*(")||2^^ < c3o}, = ^ A„ , for all * e . (54) 

?l=0 Tl = 

Taking in particular An ~ nif^^ for n G N, we obtain the number operator 

oo 

N^Y^ nl-Hr, . (55) 

n=0 
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We say that two self-adjoint operators B and C on a Hilbert space strongly commute if their 
spectral projections mutually commute. This is equivalent to the commutation of their resolvents 
for some z e C\R and also to the commutation of their associated unitary groups. More precisely, 
B and C strongly commute if and only if for all t, s G M 



Proposition A.l. Let A and N he the operators given by Il54\ )- f55\) . The following assertions 
hold: 

(i) A and N are self-adjoint. 

(ii) For any bounded Borel function on R, 



/(A) = ^/(A„) 



n=0 



(Hi) The operators A and N strongly commute. 

(iv) IfT>n is a core for An for each n G N then (BnloT^n is a core for A. 

(v) For any real polynomial p the operator A + p{N)\x>{A)nv{p{N)) is essentially self-adjoint and 

oo 

A + p{N)\-i,(A)r^V(p(N)) = X!^" • 

n=0 

Proof, (i) Clearly, ^ is a densely defined operator. It is also symmetric, since for any ^, $ e T^{^) 

oo oo 
n=0 n=0 

For any ^ £V{A) and ^ eV{A*), 

OO 

Hence the inequality holds 



n=0 



< \\A*mu II* 



By taking any -J-^") e , this means G T^iA^) = 2?(A*J = ^{An) . The extension to any 

^ G y, gives $ G 'T>{A). This proves that A and N are self-adjoint. 

(ii) For each n e N, the map t ^ e^tA^-ztA„^(n) of class for any e X'(A„) by the Stone's 
theorem with the derivative 

Hence, for any ^ G Ql?(y4„) (and then for any * G "H since o^(^n) dense in we see 
that for alH G K 

OO 

e**^* = ^e^*^"*(") . (56) 



By functional calculus we extend the identity ([56]) to any bounded Borel function / on 
(iii) By using (ii), we get for all s,t and \l/ G "H 

oo oo 
n=0 ri=0 
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(iv) The algebraic direct sum = ®°^q2?(A„) C 'D{A) is dense in H and invariant with 
respect to the group (e**'^)tgR. Therefore, Vfin is a core for A. In the other hand, the subspace 
^/in = ®n=o^n Satisfies 

Hence 2?^i„ is also a core for A since A^-pfi^ — ^\v(A)- 

(v) The operator B = X^^g^" +P(^)ll'Hn (with its natural domain) is self-adjoint by assertion 
(i). It is clear that 

Vfin = K%i^{An) c v{A) n V{p{N)) C V{B) , 

and furthermore 

B\Vfi„ = A + p{N)\T>fi„ C A + p{N)\T>{A)nv{p{N)) C B\t>{b) ■ 
Therefore, the operator A + p{N)^x){A)nv{p{N)) is essentially self-adjoint since -B|Dj,j„ = -B|d(b)- D 



B Second quantization 

For reader convenience we recall here the general framework of second quantization and put some 
related notations. 

The phase-space, a complex separable Hilbort space, is denoted by Z with the scalar product (., .). 
The symmetric Fock space over Z is defined as the following direct Hilbert sum 

oo n 

rs{z) = ^\Jz, 

n=0 

where \/" Z is the n-fold symmetric tensor product. The orthogonal projection from onto 
the closed subspace Z is given by 

■ cre6n 

Algebraic direct sums or tensor products are denoted with a "'^ superscript. Hence 

alg n 
TiGN 

denotes the subspace of vectors with a finite number of particles. The creation and annihilation 
operators a*{z) and a{z), parameterized by e > 0, are then defined by : 

They extend to closed operators and they are adjoint of one another. They also satisfy the e- 
canonical commutation relations (CCR): 

[a(zi),a*(z2)] =e(zi,Z2)ll, [a* (zi) , a* {Z2)] = = [a{zi) , a{z2)] . (57) 

The Weyl operators are given for ^; e ^ by 

W{z) = e75l''*(^)+"(^)] ^ 

and they satisfy Weyl commutation relations in the Fock space 

W{zi)W{z2) = e-^i-" W{zi + Z2), zi,Z2 e Z . (58) 
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The number operator is also parametrized by £ > 0, 

For any self-adjoint operator A : Z D — > -Z, the operator dr(A) is the self-adjoint operator 
given by 



dr(^)|^n,alg-p(^) = £ 



B.l Weyl, Anti-Wick quantized operators 

Let P denote the sot of all finite rank orthogonal projections on Z and for a given p G P let 
Lp{dz) denote the Lebesgue measure on the finite dimensional subspace pZ, with volume 1 for an 
orthonormal hypercube in pZ. A function / : Z — )• C is said cylindrical if there exists p € P and 
a function g on pZ such that /(z) = g(pz), for all z E Z. In this case we say that / is based on 
the subspace pZ. We set Scyi(Z) to be the cylindrical Schwartz space: 

(/ e Scyi{Z)) ^ (3p €F,3g€ S{pZ), f{z) = g{pz)) . 

The Fourier transform of a function / e Scyi{Z) based on the subspace pZ is defined as 

JpZ 



and its inverse Fourier transform as 



/W= / Him e2-R<=<-^«) Lp{dO. 



With any symbol h G Scyi{Z) based on pZ, a Wcyl observable can be associated according to 

bWeyi ^ f jr[j,](^) W {V2tt z) Lp{dz). (59) 

JpZ 



Notice that 6^*^^' is a well defined bounded operator on H for all b G Scyi{Z) and that this 

quantization of cylindrical symbols depends on the parameter e. 

We also recall the Anti-Wick quantization thought its usual finite dimensional relation to Weyl 
operators. In fact, we have 



LA—Wick 



' pZ 



Weyl 

' pZ (7r£/2)d™S=2 

Fibm W{^/2^,^) e-'^^^^U Lp{d^) , (61) 



for any b e S{pZ) by setting 



b * 7(^)= / b{zMz-z')Lp{dz'). 

JpZ 



B.2 Wick quantized operators 

For any p,q G N, the space Vp^q{Z) of complex- valued polynomials on Z is defined with the 
following continuity condition: 

p q 

b e Vp^giZ) iff there exists b G C{\/ Z, \/ Z) 
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such that: 

b{z) = {z®\bz®P) . 
On these spaces the norms are given by |&|p,g = \\b\\c{\/'' z,\/'' z) ■ 

The subspace of Vp q{Z) made of polynomials 6 such that & is a compact operator is denoted by 
VZi^) ■ 

The Wick monomial of a 'symbol' h £ Vp_q{Z) is the linear operator 
defined as : 



0\y^ z = Mp,+^){n) (^n-p)\ Sn-p+q[b®T[' j . (62) 

where 6(g) ll^""'') is the operator with the action (h ® = (hip'^P) Notice that 

^Wick cjgpends on the scaling parameter e . When h is an unbounded operator with domain 'D{h) 
containing Y^'"'^ T) , the formula ((62|) makes sense when applied to ^I^ e y^'^'f p . 

Proposition B.l. For b G Vp,q{2), the following number estimate holds 

(N)~H^*^'=(N)-' < |5| . (63) 

£(«) ^P'" 

An important property of our class of Wick polynomials is that a composition of b^^'^'^ ° b^^'^^ 
with 61,62 G ffipjf'^p.gl-^^) is a Wick polynomial with symbol in ®^^'Pp^q{Z) . This was checked 
with a convenient writing in [3] and widely used also in [S] [5] . 

We need some notations: For 6 £ Vp_q{Z), the fc-th differential is well defined according to 

k k 

d';b{z) £ {\J Zy and dib{z)£\JZ, 

for any fixed z £ Z . Actually (V*^ Z)* is the dual of (V*' Z) with a C-bilinear duality bracket. For 
two symbols bi £ Vp-^q^{Z), i ~ 1,2, and any fc € N, the new symbol 9^6i.9|62 is now defined by 

d% . dlb^iz) - {d%{z),d%iz))^y.^y y.^ . (64) 

We also use the following notation for multiple Poisson brackets: 

{bi,b2}^''^ = d%.dlb2- d%.d%, 

{bi,b2} = {6i,62}«. 
With these notations, the composition formula of Wick symbols has a very familiar form. 
Proposition B.2. Let 61 £ Vp^,q^{Z) and 62 £ 'Pp2.q2{Z) . 

For any k £ {0, . . . , min {pi, 52}}, 9^6i.9|62 belongs to 'Pp2-k,qi~k{Z) with the estimate 

Pi! 92! 



The formulas 



, . , , , Wick 

'min{pi,(j2} J, 

(t) bf'^^obf'^^ = I H IT 9%.dlb2 



k\ 

k=0 



(e^<^-^->6i(z)62(c.)U^, 



Wick 



f max{min{pi,g2} . min{p2,<3i}} j, 

(zi) [6f-^6f-'=] =1 J2 IT {^1=^2}^'' 



k\ 

k=l 



hold as identities on Jifm- 
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/ \ Wick 

Combined with Proposition ED and [b^^'^'^y = b{z) this also gives the 

Proposition B.3. For b G Vp,qiZ), {N)~^b^"''' and 6'^'"='^ (N) - ^ extend as bounded op- 
erators on T-L with norm smaller that Cp^q\\b\\£^(\jp z^y z) , for all e G (0,e) . 
We will also need some more particular estimates stated in the following two lemmata. 

Lemma B.4. Let A be a self-adjoint operator on Z with A > 11. For any polynomials bi e 'Pi^2{Z) 
and 62 G 7^2, 1(2) the estimates below hold true: 

ii) II (dr(A) + Vn + 6f -'^^ (dr(A) + Vn + i)-i|| < \\a-^'^ h (n ® A-^i^)\\^^y. z,z) , 
(m) ||(dr(A) + VN + i)-i bT"'' (dr(A) + \/N + i)-i|| < m®A-^'^)b2A-^/^\\^^z,y^z)- 

Remark B.5. The term \/N can be absorbed in dr(j4) + 1, if one accepts constants larger than 
1 as factors of the right-hand sides of (i) and (ii). 

Proof. The estimate (ii) follows from (i) by taking the adjoint. Let us prove (i). 
For $, * e e^'io V"''''^ we write 



n=2 

CO 

= e'^/^y/n{n-lY{{A^''^ ® jl(n-2))^(„-i) ^ ^^^jj ^ ^1/2 ^ ^(„)^ ^ 



n=2 



with 



Ba = [(A-i/26i)(ll ® ^ ]l("-2) _ 

Hence, by Cauchy-Schwarz inequality, we get 

\{^,bY'^H)\ < \\A-'/%i^[^A-'/^)\\^(y.z,z) 



1/2 



X J2 s'^Wnin-iniA'/^ ^(n-2))^(„-i)||2 



\n=2 J 

/ 00 \ 1/2 

X ( E ^'^V«(" - 1)2||(]10 €D ]l("-2))<j,(n)|j2\ _ 

Now, we observe that 

£3/2^„(„_ l)2||(Al/2 ^ ;il(n-2))^(«-l)||2 < 2 ^sl" " 1) (**""^\ (e(n " 1) ) ( A (g) ll^""^' ) ) 

< 2(^'("-i),VNdr(A)*("-i') 

and 

e'l'^^Jn{n - l)2||(ll(g) ^1/2 ^ ii(«-2))$(«)||2 < (11(g) A (g) ll*""^' 

< v^($("\ne(A® ll("-i))$(")> 

< ($("), VNdr(A)$(")) . 

On the other hand, with the inequality 2ah < + 6^, we see that 

2(5-("-i),yNdr(yi)*("-i)) < (*("-i\(N + dr(A)2)*("-i)) < ||(VN + dr(A))«'("-i)f 
and 2($("\ VNdr(A)$(")) < ($("\ (N + dr(A)2)^'(")) < || (Vn + dr(A))$(") f , 
where the last inequalities come from 2V]Vdr(A) > 0. Therefore, we obtain 

|(*,6f'^'=$)| < (11(8^-1/2)11 ||(VN + dr(A))*|| ||(VN + dr(A))$|| 

and hence the estimate extends to $, 'J e I?(v/N + dr(A)) n Hfm- This means that the operator 
(dr(A) + \/N + l)-^bY"'''{dT{A) + \/N + l)|i,\vN+dr(A))nw extends to a bounded operator 
satisfying (i). □ 
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Lemma B.6. Let A, B two self-adjoint operators on Z with VIA) C 'D{B) and B >0. Let C be 
a self-adjoint operator on Z such that 2?(C) C 'D{B2) where B2 ~B®'&+)\®B. Then the 
estimates below hold true: 

(z) ||(dr(B) + N + dT{A) {AT{B) + N + < ||(1 + B)-' A{1 + B)-'\\c(z) , 

{ii) ||(dr(B)+N2 + l)-i/2c'^-fe (dr(B)+N2 + l)-V2|| < \\{l + B2)-^'^C{l + B2)-^'^\\c(y-z) 



Proof. We follow a similar argument as in the proof of Lemma IB. 41 Indeed, Cauchy-Schwarz 
inequality gives for every 

|(vI',dr(A)$)| < \\{l + B)-^A{l + B)-^\\c(z) I ^en||((l + i?)®]l("-i))*(")|pj 

X 1/2 

xK]£n||((l + B)0]l("-i))$(")|n ■ 



Now, we observe that 



en 



< ||dr(i + s)*(")|p^„2. 



since in the sense of quadratic forms dr((l + BY) < dr(l + B)^. Hence we obtain 

|(*, dr(A) $)| < 11(1 + B)-'A (1 + B)-^ciz) ||dr(l + B)n l|dr(l + B)m . 
This proves (i). 

Expressing (7^''^'^ as a quadratic form for ^f,*!) e ©°'f g V"'"'^ £^nd then applying Cauchy- 

Schwarz inequality yield 

K^^C^-z^*)! < |l(ll + B®ll+]l®S)-i/2c(ll + B®ll+ll®S)-i/2||^^^,^^ 

/ 00 \ 1/2 

X J2 - 1)11 [(11 + S (g) 11 + 11 ® (g ]l("-2)]^(«)||2 

/ 00 \ 1/2 

X ^e2n(n-l)||[(]l + B(g)ll+]l®S)i/20]i(»-2)]$(")||2 
Due to the symmetry of the vector we remark that 

S^n^W [(11 + B (g) 11 + 11 ® (g = £2^2 ^$("), [(11 + B «) 11 + 11 (g (g ll("-2)]$(n)) 

= 2($("),(l + dr(S))N2$(")) 

< ||(i + dr(s) + N2)$(")||2 

So we obtain 

l(vi/,c^-feci>)l < ||(i + i?2)-i/2c(i + S2)-i/2||^(V^2) ||(i + dr(s) + N2)M/|| 
x|j(i + dr(B) + n2)$|| . 

This proves (ii). □ 



C Absolutely continuous curves in Prob2(2^i,M) 

This section firstly gathers results presented in [3 about Borel probability measures on a separable 
real Hilbert space which are weak solutions to continuity equations. In a second step, we shall 
adapt it to a complex Hilbert space Zi endowed with its real euclidean structure. 
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C.l Absolutely continuous curves in Prob2(-E') 

Let E he a real Hilbert space with scalar product ( , ) and norm | | . The symbol Probp(-E) (resp. 
Prob(£')) refers to the set of Borel probability measures fi on E such that \x\p dii{x) < +00 
(resp. with no momentum condition), and we simply work with p — 2. On Prob2(-E), the 2- 
Wasserstein distance is defined by 

:= min I^J \xi - X2\l: dfi{xi,X2) ; = , 

where Ilj : E^ --^ E is the natural projection, j = 1,2 . Narrow convergence of a sequence (/^„)„gN 
of Prob2(£'), with a uniform control of |a;p d^n is equivalent to the convergence on Prob2(£') 
(see Proposition 7.1.5 in [3]) . Remember also that the tightness property of subsets of Prob2(£') 
can be checked in the infinite dimensional case with the weak topology, or after introducing a 
Hilbert basis (e„)„gN', with the distance duj{xi,X2) — X^new ^^"^^T+np"^^ ■ This use of weak or 
topology, is done also when considering probability measures on set of absolutely continuous 
curves in E . 

This tightness property is called the weak tightness property in [3] since it refers to the weak 
topology on E. Especially when one considers the narrow convergence in Prob2(i?), there is 
a weak narrow convergence and a strong narrow convergence (see the discussions about this in 
Chapter 5 and 7 of [3]). The terms "narrow convergence" or "narrow continuity" refer to the 
strong ones and we shall specify "weak narrow convergence" and "weak narrow continuity" when 
necessary. 

We recall two results of 0] and give a complete proof in the infinite dimensional case of the second 
one, for the sake of completeness (it is left as an exercise to the reader in 
The following result is the second part of Theorem 8.3.1 in [3] with p = 2 . 

Proposition C.l. Let I be an open interval m M . If a -continuous curve jit ■ I ^ Prob2(i?) 
satisfies the continuity equation 

atAit + V^(wt^t) = (65) 

in the weak sense 

{dMx, t) + {vt{x) , V,(p(x, t))E) dfitix)dt = , e C^,yi{E X /) , (66) 

for some Borel velocity field Vt, with \vt\L'^(E,iit) ^ -^^(^)j then '■ I ^ Prob2(£') is absolutely 
continuous and |/i'|(i) < \vt\L^{E.fj.t) /'^'^ Lebesgue almost every t € I . Moreover for Lebesgue almost 
every t E I , Vt belongs to the closure in L'^{E, fit) of the subspace spanned by | V(/3, € Cj^^yi{E)^ . 

The previous result concerns non regular (non Lipschitz) vector fields for which there is no 
uniqueness result for the Cauchy problem. Remember that the infinite dimensional case, which re- 
lies on the cylindrical integration of vt and cylindrical disintegration of the measure /it, requires the 
introduction of such singular vector fields (see the proof of Theorem 8.3.1 in ^). Nevertheless an 
interpretation of the continuity equation (|65|) ((66|) in terms of characteristic curves can be done via 
a probabilistic representation. For the sake of completeness, we adapt the proof of Theorem 8.2.1 
stated in ^ for the finite dimensional case, to our infinite dimensional case. For T € (0,-|-cx)), 
consider the set Ft — C"([— T, T]; £') endowed with the norm |7|oo,T = maXfg[_3-^3-] |7(i)|£; . For 
a Borel probability measure rj defined on £^ x P^, consider the time dependent Borel probability 
measure /ij* defined by 

I ^dfi^l^ ( Mt)) dvix, 7) , e Cl,yi{E), t e [-T, T] . (67) 

The measure fi^ is the push-forward of r] by the evaluation map 

et : (x, 7) G £; X Ft ^ 7(i) eE, for t e [-T, T] . 
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Proposition C.2. Let /it : [— T, T] — )■ Prob2(-E) be a W'^ -continuous solution to the continuity 
equation (j65p(|66p . with I = {—T,T), for a suitable Borel vector field v{t,x) = Vt{x) such that 
\vt\L^{E,fj,t) ^ ■ Then there exists a Borel probability measure r] in E x Tt such that 

(i) r) is concentrated on the set of pairs (a;, 7) such that 7 G AC'^{[—T,T];E) is a solution to the 

ODE 7(i) = vti^it)) for Lebesgue almost every t G (— T, T) with 7(0) — x ; 

(ii) /it — for any t G [— T, T], wif/i /ij* defined as in ([67]). 
Conversely, any r) satisfying (i) and 




\vt{l{t))\E dr]{x,^)dt < +00 , 



induces via (j67p a solution to the continuity equation, with /io = 7(0)^77 . 

Remark C.3. The notation AC'^ ([~T , T]; E) refers to the set of absolutely continuous curves 
in E with L^([— T, T]; _E) derivative. We keep the notation $,t/i of differential geometry, for the 
push-forward or direct image of a measure pL, by the Borel map $ . 

Proof. The result is proved in [3] wtien E is finite dimensional. The proof of the second (converse) 
part of the statement, is exactly the same as in finite dimension, after replacing regular (Lipschitz) 
test functions by cylindrical ones. We now show, for the first part, how the infinite dimensional case 
is deduced from the finite dimensional result, following an approximation scheme like in the proof 
of [31 Theorem 8.3.1]. After introducing an Hilbert basis (e„)„gN* of E, the maps ir^^ : E ^ W^, 
^d,T : i^ti _^ £; and tt"^ : E ^ E arc defined according to 

w'^ix) = {{ei,x),...,{ed,x)), 

d 

7r'^''^(yi,...,yd) = ^Vje-j, 

;;^d d,T ^ d 

With the measure fit G Prob2(i?), the measure fif G Prob2(M'^) is defined by fif = Tr^/ij and 
{f^t.y, y denotes the disintegration of /it w.r.t fif . Within the space E endowed with the basis 

(e„)„gN* , fit is nothing but (g) Sq in the decomposition Z = F^x with Fd — span(ei, . . . , e^) . 
The vector field vf (resp. vf) is defined on (resp. on E) by 

vf{y)= [ TT^ix) dfit.yix) , yeR'' 

resp. vf{y)^ / %'^vt(x) dpLt^y{x) , yeE. 

J(7r<i)-i(Tr''iy) 

Within the proof of Theorem 8.3.1 in it was checked that (resp. fif) is a weak solution to 
the continuity equation 

9tMf + V^(t;fMf) = 0, 
resp. 5tAf + V^(f)fAf) = 0. 

with the following properties: 

1) \vt\L^(E,fif) = L2(Rd,^d) < \\vt\\L^{E,dtit) ; 

2) W'ifii,fii) < III \<\L-{^-,i^t) < It'" h\L^(E,^,) dt, for -T<t,<t2<T- 

3) the sequence {fif)d£W converges weakly narrowly to /it with the estimate 

W^{fit,,tit,) < liminfiy2(^rf dj < f ' \y^\^^ dt, -T<ti<t2< T. 
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The set of continuous maps from [-T, T] to R"^ is denoted by . By using the finite dimensional 
result, stated in Theorem 8.2.1 of there exists for any d € N a probability measure, rj'^ , on 
R'' X such that the properties (i) and (ii) hold when {iJ,t,Vt,E) is replaced by {^f , vf ,9."^) . 
Equivalently the result can be formulated in E after using {ilf,vf,E) instead of (/if , , K'') and 
using fj'^ = {it 
which satisfy 



i^'^''^)*'nt ■ Hence we have a sequence {ri'^)de'M of probability measures on ii^ x Ft 



^ o ^'^(7(i)) df,\x, 7) , V(p e C°(M'*), t e [-T, T] 



ExFt 



where ip o tt'^ can be replaced \yy <f o ■k'^ with ip e C^{Fd) ■ 

After some regularization done |3| ppl79-180, it is proved that any measure r)** satisfies 



|7(0I' dtdi)'^ < 



T 



T 
'T JE 



vt{x)\'^ dfitdt. 



(68) 



Since the functional g J^j, \g(t)\'^ dt , defined on {g € Tt , g{0) — 0} and set to +00 if 5 ^ 
AC'^{[—T, T];E), has compact sublevel sets in Ft , the two mappings 

: {x,-f) e E X Tt ^ X e E , r'^ : (x,^) e E x Tt ^ g^,x = 7 - x € Ft 

give rise to (weakly) tight families of marginals {rlr)'^) den = (Ao)deN and {r'^r]'^)den ■ Remember 
that the compactness of subsets of E or Ft is considered with the weak topology on E or the 
distance . Hence the family (r)'^)^^^ is (weakly) tight in Prob(i? x Ft) and we take for rj a weak 
narrow limit point of i)'^ . By assuming the test function ip in (|68p to depend only on d' coordinates 
with d' < d, and by taking the limit d — > +00 while d' and p are fixed, we get 



{ipon'^ ) dfit 



((po7r'^')(7(i)) dr]{x,-f), 



Exr 



for all p> e C^i^'^') and t e [-T, T], where (p ott'^ can then be replaced by any cylindrical function 
or Borel bounded function on £^ . It remains to prove the condition (i) for rj, namely that this 
measure is concentrated on curves verifying 7 = v{^{t)) for Lebesgue almost every t g {—T,T) 
(7(0) = x is already known) . 

The estimate (8.2.6) used in [3l for the finite dimensional case, provides the inequality 



Ws{-f{s)) ds\^ drj'*(x,7) < {2T) 



-T 



\vf-wt\'d^if dt, 



for any family Ws{x) = 10(8, x) uniformly bounded continuous function from [—T,T] 



to 



After assuming that w actually belongs to C^([— T, T] 



pd' . Tad' 



) with d' < d fixed and, by using 



Wt = n 



d'T , 



Wt o TTd' e C"([— T, T] X E]E) , taking the limit as d — )■ 00 gives 



Wshis)) dsp dT]{x,j) < (2r)limsup 

d^oo 



I ExTt 

But the condition 1) for \vf \ 1^2 1^^,(1. f) is easily extended to 

\vt - Wt\mE.p.f) < \'Vt - Wt\L^(^E,t,t) 

by the same argument, relying on 



|5f -tutpdAf dt. 



(■of - Wt , x) dfi 



< 


L 






< 


vt - 



{%d{vt{x) - Wt{x)) , x{t^'^{x))) dutix) 



^X&L\Er^t). 
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This uniform upper bound leads to 



/ \-i{t)-x- ( Ws{j{s)) ds\^ drj{x,-f) <i2T) f f \vt{x) ~ Wt{^)\^ dfitdt . 
JexFt Jo Jo Je 

According to the last statement of Proposition IC.ll vt can be approximated in L^{E,^t) by a 
sequence of bounded regular cylindrical functions, (i()t,n)ngN ■ By possibly truncating with respect 
to times t — >■ Wt^n so that \vt — Wt,n\L^(E,fit) — 1 ^r a.e. t £ {—T,T) and all n S N*, Lebesgue's 
Theorem implies 

hit)~x~ f «.(7(s)) ds\^ drj\x,^)=0, 
exTt Jo 

which ends the proof. □ 
Below is a consequence of the above probabilistic interpretation when the Cauchy problem jit) = 
''^tilit)), 7(0) = X admits a unique solution for all x € E . The fact that we have to pass by the 
probabilistic representation is a real question. Contrary to the finite dimensional case, the well- 
posedness of the Cauchy problem, even with the standard Picard's contraction argument, defining 
a flow on the whole space E, does not give a representation formula for observables. The point is 
that the natural observables, or test functions, are cylindrical functions, a property which is not 
preserved a nonlinear flow. 

Proposition C.4. Let /it : R — > Prob2(-E) be a -continuous solution to the continuity equation 
(|65p(|66p for a suitable Borel velocity field v{t,x) = Vt{x) such that \vt\L^(^E.fn) ^ ^^([^^i^]) 
all T > . Assume additionally that the Cauchy problem 

at) = M"/{t)) , jis) = x 

or ^ 

7(t) = ^ + '"silis)) ds, 

admits a unique global continuous solution on M for all s Cz R and all x £ E, such that j(t) ~ 
<i>(t, 5)7(3) defines a Borel flow on E (i.e. <l>(i, s) : E ^ E is a Borel function for all t, s £ 
Then the measure fit satisfies 

V<,seR, /it = $(i,s),/i,, . 

Proof. It suffices to work with t £ [—T, T] as in Proposition IC.2I Since the evaluation map 
Ct : E X (a;,7)rT 7(i) £ E is a continuous, thus Borel, map. The relation fit = fJ-t defined 
according to (|67p extends to any bounded Borel function (p on E: 

ip dfit= / v{i{t)) dvixn) ■ 

JExFt 

By using 7(t) ~ <i>(t, 5)7(5), with <i>(t, s) Borel, we deduce 



ipdfit^ / [^°'^it,s)]{"f{.s)) dr]{x,"f) ^ / [if o ^(t, s)] dfi, 
J exVt Je 

which is nothing but fit = <i>(t, s)^,fis . □ 



C.2 Application to Hamiltonian fields 

We finally specify how these results apply to our case, when the phase-space Zq is a complex 
Hilbert space and the velocity field is associated with a (singular) Hamiltonian vector field, only 
defined on Zi C Zq . 

Consider a complex Hilbert triple Zi d Zq £1 Z^i, with Zi densely continuously embedded in Zq 
and -E_i being the dual of Zi for the duality bracket extending (zi , 22)2:0 ■ The dual of a complex 
Hilbert space Z while keeping the C-bilinear duality bracket, written u .v in ([M]) . is still denoted 
by Z* . In the case treated in the article Zq = L'^{R'^, dx), Zi = ^^(R'*) and Z_i = i?-i(M'^) . The 
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space Zq is endowed with its scalar product {zi , 2^2)20; ^^^^ euclidean structure with Re {zi , Z2)zo 
and its symplectic structure a{zi,Z2) = Im {zi , Z2)zo ■ On Zi we will use the hermitian {zi , Z2}zi 
and euclidean scalar product 

(-21 , Z2)zi,R = Re {Zi , Z2)zi ■ 

For a cylindrical function / e Scyi{Zo), based on pZo, the differentials dzf{z) and dzf{z) are 
defined 

J pZo 

Hence dzf{z) is a continuous C-linear form on Zo while dzf{z) € ^0 • This notation is coherent 
with the definition of dzb{z) and dzb{z) when 6 is a Wick symbol in ®p''^q'Pp^q{Z) . 
A function / S Scyi{Zi) is given by 

f{z) = ip{{^i , z}z^,- . . , (€jv , z}z^) = , z},...,{r]N , z)) 

with . . . wm) € 5(R^) and ^1, . . . , € Zi and . . . ,7]^ € such that {^j , z)zi = 

{r]j , z) for all z & Zx . The derivatives d^f and d^f are thus given by 

N 

'iz&Zi, dzf{z) = ^dw.(p{{f]i, z),...,{r]N, z)){r]j\ e Zl^ , 

N 

yz&Zi, dzf{z) = ^dyj-ip{{r]i, z),...,{t]n, z))\r]j) G Z_i . 

When h{z) is an unbounded polynomial on Zq but which happens to be a real-valued Frechet 
C^-function on Zi, the derivatives dzh{z) and dzh{z) are defined only for z G Zi and we have 

yzGZi, dzh{z) e Zx , dzh{z) e ^* i . 

When / e >Scy;(2o) (resp. g G 5cyi(2i) or h) is real valued differentiating f{z + te) a.t t = 0, t £ M, 
for any e G Zq (resp. any e e ^1) leads to 

dzf{z).u = {u , dzf{z)) , zeZo,ueZo (69) 
dz9{z).u = {u , dzg{z)) , z e Zi , u G Zi . (70) 
dzh{z).u = {u , dzh{z)) , z G Zi,u G Z-\ . (71) 

Note that the Poisson bracket 

i {h, b} (z) = i {dzh.dzb — dzb.dzh) (z) , z G Z\ 

is well defined for h G Scyi iZ\ ; M) and our aim is to write it as the real scalar product 

(t;(^), (V6)(^))2i,M, zgZ^. 

Definition C.5. For a cylindrical function on Zi, f G Scyi{Z{), the gradients and V are 
defined by 

yzGZ,,uGZu {u, Vzf{z))z, = {u , dzf{z)) , 
V = 2Vj. 

Remark C.6. • Although it is not necessary, these definitions can be justified by introducing a 
complex conjugation u ^ u on Zp, which re mains a conjugation on Z\ , that is an isometric 
C-antilinear application such that {u,v)zo,i = {u, v)zo,i ■ When Zq = L'^{R'^;C) and Z\ 



H^{M.'^) this is the usual pointwise complex conjugation. 
For real valued function one then sets 

Vij/ = V^/ + V,/ and Wif = \ (V^/ - V,/) . 



so that 



V^f = ^{VRf + iVif) . 



Xr 
Xi 



Similarly, an element X ofZ\ can he written X = Xji+iXj withXjij = Xrj orX = 
and the real scalar product 

{X , Y)z,,M = Re{X, Y)z, = {Xr , Yr)z, + {Xi , Yi)z, • 

Then the definition of the gradient of a real cylindrical function f becomes 

• It is important to notice that we do not use the Zi-gradient for the real valued function 
h{z) , hut keep the derivative, dzh{z) modeled on the duality hracket { , ) ■ With a complex 
conjugation and since h is real valued, it can he decomposed into dzh = 5 {dRh + idih) and 

1 i 
—idzh = -djh dph . 

2 2 

Lemma C.7. With the ahove notations and assumptions the equality 

yz€Zi, i {h, b} (z) = 2Re {-id^h{z) , V^b{z))z, = {v{z) , Vb{z))z,,M , 
holds for any b G Scyi{Zi;R) with v{z) = —idzh{z) . 
Proof. It suffices to compute 

i{h,b} = i[dzh.dzb — dzb.dzh] 



{dzb , dzh)zo - {dzh , dzb)zo 
-2Im {dzh, dzb) = 2Re {-idzh, dzb) = 2Re {-idzh, Vzb)zi 
{-idzh, Vb)zi,v. ■ 



□ 



Proposition C.8. Let Zi c Zq c Z-i he a Hilbert triple of separable complex Hilbert spaces. 

Consider a time dependent real sesquilinear form, z — > h(z,t) on Z\ which is Frechet-C^ and such 
that Zi xRb {z,t) — >■ {dzh{z,t) , dzh{z,t)) G Zi x Z^-^ is strongly continuous. Assume also that 
the time-dependent Hamilton equation 

idtzt = dzh{zt, zt, t) , zt=s = z 

admits a unique continuous solution Zt = <^{t,s)z for all s € M and all z & Zi, with ^{t,s) : 
Z\ — >■ Z\ Borel. 

Consider a time dependent measure e Prob2(^i) which satisfies 

• t ^ G Prob2(^i) is -continuous. 

• For allT> 0, \dzh{t)\L2^z,f.,) € LH[-T,T]) . 
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• The time- dependent probability measure fXt is a weak solution to 

dtfi + i {h{t) , ^} = , 
namely for all ip € Cf^cyii^^ ^ R), 

{dt(p{z, t) + i {h, (fi} {z, t)) dfit(z)dt = . 



Then the measure fj,t satisfies 

Vt, seR, ^it ^ <^{t,s)*fis 

and it is unique when /iq *s fixed. 

Proof. We apply Proposition IC .41 while E — Zi is endowed with its euclidean structure (zi , Z2)ziM = 
Rc {zi , Z2)zi ■ Lemma I C. 71 savs that the weak Liouville equation is 



yipec^,yi{z,xR), / / idMz,t) + {v,Vip)it,z)) d^it{z)dt^o, 

Jm Jzi 

with v{z,t) — —idzh{z,t) . The measure fit is a weak solution to 

dtfji + ^^{vfi) = 

where V and V"^ are defined according to the real structure on Zi . Our hypothesis on fi and h 
cover all the assumptions of Proposition IC.4I □ 



D Weak conditions for the potential V 

Let L^d be the Lebesgue measure on R''. Let < p < cx), a Lebesgue measurable function 
/ : R'' — )■ C is said to be weak integrable, or shortly in L'P'°°{W^), if there exists a constant 
c > such that for alH > 

L^.{x:\f{x)\>t}<d-/t^ . 
Two functions in U''°°{K^) are equal if they are equal Lgd-almost everywhere. The quantity 

ll/L.oo = inf{c:LR4:E: >^}<cV^^V^>0} 

= sMtL^.{x : > tY'^} 

t>o 

defines a complete quasi-norm on LP'°°(R'') with ||/||p,oo < ll/llp- 

By combining Hunt and Marcinkiewicz interpolation theorems according to [34| , l55 [ fTl ] ). the Young 
and Holder inequalities can be extended to weak spaces. 

Proposition D.l (generalized Young's inequality). Let 1 < p,q,r < oo such that ^ + ^ = 1 + f • 
There exists a constant Cp,q > such that for all f G ^^(R'') and g G L'^'°°(R'*) 



||/*5l|r < Cp,g ll/llp ||g|| 



q,oo • 



Proposition D.2 (generalized Holder inequality). Let 1 < p, q,r < oo satisfying ^ + ^ = 7- There 
exists a constant Cp^q such that for all f G LP'°°{M.'^) and g G L'(R'*) 

< Cp,q ||/||p,oo hWq- 

Proposition D.3 (Hardy inequality). Suppose that d > 3 and V G L'^'°°{M.'^). There exists a 
constant c > such that for all u G H^(M.'^) 

\\Vu\\2 < c||u||^l(Rd) . 



44 



Proof. For u e L2(M'^), we can write (1 - A)-i/2u(a;) = G * u{x) with G the inverse Fourier 
transform of (1 + It is not difficuh to prove that G e L^'°° (see [S51 Exercice 50]). 

Hence, we conclude that 

, Holder 

\\Vil-A)-'/^uh^\\VG^uh < Gi\\V\\d,^\\G^u\\^ 

Young 

< G2\\V\Uoo\\G\\^,^\\uh 

□ 

The above proposition provides a class of potentials which are bounded multiplication operators 
from i?^(R'*) into L^(R) when the dimension d > 3. For lower dimension, the Sobolev embeddings 
give at once: 

• if d= 1, y e L'^{R) + L'^iR) then V G C{H\R) , L'^ (R)) . 

• if d= 2, y e LP(M2) + for p > 2, then V £ CiH^iR"^), L^iR"^)). 

We denote by LP{R'^) + L^{R'^) the space of Lebesgue measurable functions / such that there 
exists (/„)„eN G LP{R'^ f satisfying lim„_>oo ||/ - /n||oo = 0. 

Lemma D.4. For < p < q, 

L«^°°(M'^) c LP{R'^) + L^{R'^) 

Proof. For e > 0, we decompose each / £ L'^'°°(R'^) into a sum f = fc + such that ~ /l|/|>e 
and = /l|/|<e- We observe that for any e > 



= p / tP-'L^.{x:\Mx)\>t}dt 
Jo 

= p tP-'L^4x:\f{x)\>t}dt + ePLfA^:\f{x)\>e} 



< cj — dt + ePLiit4x ■.\f{x)\> e} < CO 

Moreover, when e — >■ 

ll/ico = ||/l|/|<Jco <e^0. 
Therefore, each / e L'^^°°{R'^) belongs to the space LP{R'^) + L^{R'^). □ 

Proposition D.5. For any V E L^iR'^) + Lff (R'') such that V{1 - A)-!/^ g C{L^{R'^)) the 
operator {1 — A)^^^^V{1 —■ A)^^/^ is compact (i.e.: V is a relatively compact perturbation of —A). 

Proof. Let g{^) = (1 + and gm{0 — ^[o,m]i\i\)9iO- The following norm convergence 

holds 

lim grn{D)Vgm{D) = g{D) V g{D) , 

using the fact that lim^^oo \\9rn{D) - .g(-D)||£(L2(Rd)) = and .9(-D)|l£(L2(E<i)) < oo. 

By Lemma [0.41 there exist Vn E L^(M'^) such that lim„_5.oo ||^ ^ Killoo = 0. We observe now that 

the Hilbcrt-Schmidt norm of gjn{D)Vn{x)gm{D) is 

\\9niiD)Vngm{D)\\c2(L^(Rd)) = llffmlll < OO. 

Therefore by norm convergence, the operator g{D) V g{D) is compact. □ 
Corollary D.6. The potential V satisfies the assumptions in the following cases: 



ifd=l andV E L-^(R) +iff (M) 
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• ifd = 2andVe LP{M.^) + Lg°(R2) with p>2, 

• ifd>3andVe L'^^°^{R'^). 

Proof. Combine Proposition ID.31 Lemma ID. 41 and Proposition ID.51 with the observation in 
dimension d = 2 

LPiR^) + L^{R^) C L2(M2) + Ljf (M^) for p > 2 . 

□ 

In particular, in dimension d = 3 the Coulomb potential V{x) = ijjj satisfies the assumptions 
(42), (411) and (421) because A e L^^°°iR^) . 
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